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We consider random matrices of the form H = W + W , where is a Wigner matrix and V a random 
diagonal matrix. We assume subexponential decay for the matrix entries of W and we choose A ~ 1 so that 
the eigenvalues of W and XV are of the same order in the bulk of the spectrum. In this paper, we prove for a 
large class of diagonal matrices V that the local deformed semicircle law holds for H, which is an analogous 
' result to the local semicircle law for Wigner matrices. We also prove complete delocalization of eigenvectors 



and other results about the positions of eigenvalues. 
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§ ! 1. Introduction 

cn 

Consider a large matrix whose entries are random variables. On one extreme are standard Wigner matrices 
^ ■ with i.i.d. random entries, whose eigenvalues are highly correlated. It is widely believed that the behaviour of 
the eigenvalue correlation functions exhibits universality, i.e., under suitable conditions, the microscopic statistics 
depends only on the symmetry class of the ensemble. Recently, the Wigner matrices have been studied intensively 
d . and many important results have been proven [15, 16, 17, 18, 19, 21, 22, 23, 12, 13, 20, 37, 38]. 

Poisson statistics for systems represents the other extreme. It corresponds to diagonal matrices with i.i.d. 
random entries. While the eigenvalues of the Wigner matrix are strongly correlated, the diagonal randomness 
makes eigenvalues independent, hence uncorrelated. Physically, the diagonal matrix may represent an on-site 
random potential on a lattice system. Compared to the mean-field nature of the Wigner matrix, which is in 
the weak disorder- or the delocalization regime, the diagonal randomness also provides a good example in the 
strong disorder- or the localization regime. It is conjectured that, after quantization, classically integrable systems 
correspond to Poisson statistics whereas classically chaotic systems correspond to random matrix statistics. In 
terms of quantum chaos, the diagonal matrix describes the 'regular' part, while the Wigner matrix is a good 
model for the 'chaotic' part. 
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It is thus natural to consider the interpolation of the two, i.e., the N x N random matrix 

H = \V + W, (1.1) 

where V is a diagonal random matrix, or a 'random potential', and W a standard Wigner matrix. Here, W is 
properly normalized so that the typical eigenvalues of V and W are of the same order. (See Definition 2.1 for the 
precise definition of H.) The parameter A determines the relative strength of each part in this model. 

For A ^ 1 the eigenvalue density is not solely determined by ^ or 14^ in the limit TV — > oo, but can be 
described by a functional equation for the Stieltjes transforms of the limiting eigenvalue distributions of V and W; 
see [33, 34]. In general, this eigenvalue distribution, referred to as the deformed semicircle law, is different from the 
semicircle distribution. The equal strength of V and W makes it non-trivial to find the nature of the interpolation 
H. For example, the eigenvectors are completely delocalized for W whereas they are localized for V, hence the 
eigenvector localization/delocalization problem requires deep investigation of the model. We remark that there 
are some results related to this model [8, 28, 4]. 

When W belongs to the Gaussian Unitary Ensemble (GUE), H is called the deformed GUE, and it can 
describe Dyson Brownian motion [10] on the real line; see, e.g., [26]. There have been many important works 
with various scales of A: Related to symmetry-breaking, transition statistics for eigenvalues in the bulk, especially 
the nearest neighbor spacing, was studied in [32, 24] for A ~ N^^^. In this situation, the diagonal part XV controls 
the average density, while the GUE part induces fluctuation of eigenvalues. For A < 1, it was shown in [35] that 
universality of eigenvalue correlation functions holds in the bulk of the spectrum. Concerning the edge behaviour, 
it was shown in [27] that the transition from the Tracy- Widom to the standard Gaussian distribution occurs on 
the scale A N^^^^. For A <^ N~^/^, the Tracy- Widom distribution for the edge eigenvalues was established 
in [36]. To our knowledge, these results have not been established for the deformed GOE- or Wigner ensembles, 
but are expected to hold true for these ensembles as well. 

In this paper, we prove that many important properties of the Wigner ensemble, including local deformed 
semicircle law, dclocalization of eigenvectors and rigidity of eigenvalues, remain valid also for the random matrix H 
in (1.1) with appropriate modifications. The main difficulty lies in i. the control of a fluctuation originated from 
the diagonal part V; and ii. the possible difference between the semicircle law and the deformed semicircle law 
at the edge. In order to overcome the fluctuation problem, we define a random variable that only depends on the 
diagonal part and provides a leading correction term for the difference between the Stieltjes transforms of the 
empirical density and the deformed semicircle law; see Theorem 2.11. The proof of the local deformed semicircle 
law also requires the square-root type behaviour of the deformed semicircle law, which forces interesting conditions 
on V and A; see Lemma 2.6 and the Appendix. 

The paper is organized as follows: In Section 2, we introduce the precise definition of the model and state 
the main results. In Section 3, the weak deformed semicircle law is proven, which provides an initial estimate 
to prove the strong deformed semicircle law, as well as a proof for the eigenvector dclocalization. In Section 4, 
with the aid of the weak deformed semicircle law and the fluctuation lemma (Lemma 4.1), we prove the strong 
deformed semicircle law. In Section 5, we identify the leading correction term for the fluctuation. In Section 6, we 
establishes estimates on the density of states and the rigidity of eigenvalues. Technical details about the square 
root behaviour and the stability bounds for the deformed semicircle law are given in the Appendix. 

Acknowledgements: We thank Horng-Tzcr Yau for suggesting this problem to us and for numerous helpful 
discussions. We are grateful to Antti Knowles for discussions and many useful remarks. We also thank Alex 
Bloemendal, Paul Bourgade, and Laszlo Erdos for helpful comments. 



2. Definition and Results 

In this section, we define our model and state our main results. 

2.1. Free convolution. As first shown in [33, 34] the limiting spectral distribution of the interpolating model (1.1) 
is given by the (additive) free convolution measure of ^, the limiting distribution of the entries of V, and jigc, the 
semicircular measure of radius 2. In a more general setting, the free convolution measure, fii ffl /i2, of two prob- 
ability measures /ii and /i2, is defined as the distribution of the sum of two freely independent non-commutative 
random variables, having distributions /ii, /i2 respectively; we refer to [40, 30, 25, 1]. The (additive) free convo- 
lution may also be described in terms of the Stieltjes transform: Let /x be a probability measure on R, then we 
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define the Stieltjes transform of /z by 

m,{z):=[^, zeC+. (2.1) 

Note that mn{z) is an analytic function in the upper half plane, satisfying limy_).oo iy'm^{iy) = 1. As first shown 
in [39, 6], the free convolution has the following property: Denote by m^^, m^j, m^^H^^' Stieltjes transforms 
of Hi, /X2, /xi ffl /X2, respectively. Then there exist two analytic functions Wi, LO2, from C+ to C+, satisfying 
limj,_).(X) uji{iy)/iy = 1, (i = 1,2), such that 

m^,^s^l2i^) = m^i{i^i{z)) = m^,(a;2(^)) , 
LJi{z) + UJ2{Z) = z --, (2.2) 

for z e C"*". The functions oji are referred to as subordination functions. Note that (2.2) also shows that 
/Ui ffl jU2 = M2 ffl Ml • It was pointed out in [9, 3] that the system (2.2) may be used as an alternative definition of 
the free convolution. In particular, given /ii, /i2, the system (2.2) has a imique solution (m^^H^^^j "^i; ■^2)- 

The system (2.2) has been used in [34] to exploit the limiting eigenvalue distributions for random matrices of 
the form A + UBU* , with A, B deterministic or random N x N matrices and U a, N x N random Haar unitary 
matrix. Free probability theory turned out to be a natural setting for studying global laws for such ensembles; 
see, e.g., [40, 1]. For more recent treatments, including local laws, we refer to [28, 8, 4]. 

In case we choose the measure ^2 as the standard semicircular law Hsc{<iE) = ^\/4 — E'^l^_2^2]{E)dE. A 
simple computation reveals that the Stieltjes transform m^^^ = nisc satisfies 

msc{z) = , . . , z gC+ . 

z + msc{z) 

Using this information, we can reduce the system (2.2), to the self-consistent equation 

mfciz) = I , ^ g C+ , (2.3) 

J X- z-mfc{z) 

with limj,_>.co W'nT'fciiy) — 1, where we have abbreviated [i = /ii. Equation (2.3) is often called the Pastur relation. 
A slightly modified version of the functional equation (2.3) is the starting point of the analysis in [33] and also of 
the present paper; see (2.6). 

The (unique) solution of (2.3) has first been studied in details in [5]. In particular, it has been shown that 
limsup^N^Q I Imm/c(£' + i??)] < 00, E € M., and hence the free convolution measure /x/c = /x ffl j^sc is absolutely 
continuous (for simplicity we denote the density also with fijc) and we conclude from the Stieltjes inversion 
formula that 

lj,fc{E) = lim 1mm fc{E + 17]) , EgR. 

Moreover, it was shown in [5] that the density /x/c is analytic in the interior of the support of /x/c. We refer to, 
e.g., [2] for further results on the regularity of the free convolution. 

2.2. Assumptions. In this section, we define the model (1.1) in details and list our main assumptions. 
2.2.1. Definition of the model. 

Definition 2.1. Let W be an N x N random matrix, whose entries, (wij), are independent, up to the symmetry 
constraint Wij =Wji, centered, complex random variables with variance and subexponential decay, i.e., 

P (ViVjwy I > < Coe-^'^' , (2.4) 

for some positive constants Co and 6 > 1. In particular. 
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Let V = {vi) he an N X N diagonal random matrix, whose entries {vi) are real, centered, i.i.d. random 
variables with law fi. More assumptions on will be stated below. For A e M we consider the random matrix 

H = {h,j) := XV + W . (2.5) 

In the next sections, we will choose /x, such that supp/;^ = [—1, 1], hence \vi\ < 1 (almost surely), but we observe 
that varying A is equivalent to changing the support of /i. 

We define the resolvent, or Green function, G{z), and the averaged Green function, m{z), of H by 

G{z) = {G,j{z)) := ^y^\y_^ , rn{z) := ^ TrG(^) , z G C+ . 

Frequently, wc abbreviate G = G{z), m = m(z), etc. 

2.2.2. Free convolution. Following the discussion in Subsection 2.1, we define m^^ as the solution to 




with limy_>.oo i2/TOjc(iy) = 1 (or, equivalently, Immj^{z) > 0, z £ C"*"). We denote by fXj^ the corresponding 
probability measure. For simplicity, we discard the superscript A from our notation. Let us list some easy 
examples: 

i. Choosing ij, = 6i, one directly sees that fife is a semicircle law of radius 2 centered at A. 

ii. For the choice ji = ^{S-i + 5i), (2.6) reduces to a cubic equation and the support of the measure /i/c can 
be inferred from a simple analysis of the discriminant of that equation. As it turns out, the support of /i 
consists of a single interval for A < 1 and of two intervals for A > 1. For simplicity, we will exclude the 
possibility of fc having support on several disjoint intervals in the following. However, some of our results 
can be generalized to this setting. 

Hi. If /i is the standard Gaussian measure, no closed expression for /iyc exists, but the moments of fife can 
be computed recursively; see [7]. Moreover, the density of /x/c is a smooth function with Gaussian tails. 
Although the Gaussian case is important, we will not deal with measures of unbounded support, but 
comment on the Gaussian case in Remarks 2.10 and 2.14. 

2.2.3. Assumptions on XV and ji. We state our assumptions on n the distributions of the entries {vi) of V. From 
now on, we choose ^ with fx = supp[— 1, 1]. Depending on the size of the 'perturbation' parameter A, we have to 

distinguish two cases: If |A| < 1, we will assume the following: 

Assumption 2.2. [Small A] The entries of the diagonal matrix V = {vi), are bounded, centered i.i.d. random 
variables. For |A| < 1, we assume that the distribution of (uj) has a continuous density n{v), such that ^{v) > 0, 
V e (-1, 1) and n{v) =0,v^ [-1, 1]. 

This assumption ensures that the deformed semicircle law fife is supported on a single interval [Li,L2]5 with 
a square root behaviour at the edges. More precisely, wc have the following result: 

Lemma 2.3. Let A < 1 and assume that jj satisfies Assumption 2.2. Then there are —oo < Li < < L2 < 00, 
such that supp /i/c = [Li,L2\. Moreover, denoting by ke the distance to the endpoints of the support of fife, i.e., 

HE ■■= ^i^{\E - Li\,\E - L2\} , ^eK, (2.7) 

there exists C > 1 such that 

C-^^<fi{E)<C^, E&[Li,L2\. (2.8) 

In a slightly different setting this lemma has been proven in [36], see also [5, 31]. In the Appendix we explain 
how to adopt the proof in [36] to our setting. 

Remark 2.4. Above, we have chosen A independent of N . However, we may choose A = CN^^ , for some positive 
constants C and S. In this setting all our results hold true as well, in particular, in all the bounds one may simply 
replace A by CN'^. 
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Remark 2.5. Determining the endpoints Li and L2 of the support of /i/c expUcitly is, in general, not pos- 
sible, since it involves solving an implicit equation. However, for A sufficiently small, one can show that 

Li = -2V1 + A2 + 0(\^) and L2 2Vl + A2 + ^(A^). Also the measure /i/c is C'(A)-closc to the the semicir- 
cular measure of radius 2\/l + A^ in an appropriate distance, but we refrain from going into the details of this 
'perturbative approach'. 

For |A| > 1, we have to strengthen the above assumptions, since the square root behaviour at the endpoint 
of the support of may fail, for A large enough. We call a probability measure, /U, Jacobi, if it is given by a 
density of the form 

fi{v) = Z-\l + vni - , (2.9) 

where d G C^([— 1, 1]), with d{v) > 0, v G [—1, 1], — 1 < a, b < 00, and Z a normalization constant. 
We have the following result. 

Lemma 2.6. Let ^ be a Jacobi measure; see (2.9). Then, for any A e R, there are —00 < ii < < 1^2 < 00, 
such that supp fifc = [^1,^2]- Moreover, 

i. for —1 < a, b < 1, for any A G R, /i/c has the square root behaviour (2.8); 
ii. for 1 < a, b < 00, there exists Ai = Ai(/^) > 1 and A2 = A2(/i) > 1 such that 

iia. for |A| < Ai, |A| < A2. /i/c has the square root behaviour at both endpoints; 

iib. for |A| < Ai, |A| > A2, /x/c has the square root behaviour at the lower endpoint of the support (i.e., for 
E e [Z/1,0]), but there is C >1, such that 

C-\L2 - E)^ < HfciE) < C{L2 - E)^ , Ee[0,L2]. (2.10) 

Analogue statements hold for |A| > Ai, |A| < A2, etc. 

The proof of the lemma is given in the Appendix. 

For our methods to work, we have to exclude situation iib of Lemma 2.6. For |A| > 1, we will thus assume: 

Assumption 2.7. [Large A] The entries of the diagonal matrix V = (vi), arc bounded, centered i.i.d. random 
variables. For |A| > 1, we assume that the distribution of the {vi) is given by a Jacobi measure, and A and a, b 
are chosen as in i or iia of Lemma 2.6. 

2.2.4. Notations and Conventions. To state our main results, we need some more notations and conventions. For 
high probability estimates we use two parameters ^ = ^jv and ip = (Pn: We assume that 

ao <^ <^ologlogAr, ^={logNf, (2.11) 

for some fixed constants ao > 2, ^0 > 10, C > 1. These constants are chosen such that large deviation estimates 
in Lemma 3.5 hold. They only depend on and Cq in (2.4) and will be kept fixed in the following. 

Definition 2.8. We say an event f2 has {S^,v)-high probability, if 

for N sufficiently large. 

Similarly, for a given event fio we say an event fl holds with {S^,v)-high probability on Q.^, if 

P(17onn'=) <e-''(i°g^)% 

for N sufficiently large. 

For brevity, we occasionally say an event holds with high probability, when we mean (^, :/)-high probability. 
We do not keep track of the explicit value of v in the following, allowing v to decrease from line to line such that 
V > Q. From our proof it becomes apparent that such reductions occur only finitely many times. 

We use the symbols 0{-) and o( • ) for the standard big-0 and little-o notation. The notations O , o, <C 
always refer to the limit N 00. Here a <^ b means a = o{b). We use c and C to denote positive constants 
that do not depend on N, usually with the convention c < C. Their value may change from line to line. Finally, 
we write a ~ &, if there is C > 1 such that C~^|&| < \a\ < C\b\, and, occasionally, we write for A^-dependent 
quantities ajv ^ bN, if there exist constants C, c > such that |ajv| < C(i^jv)'^^|&Ar|- 
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2.3. Results. In this subsection we state our main results. 

Since we choose the measure /U to be centered, we may assume that A > 0, without loss of generality in the 
following. Fix some Aq > 0, then we assume that the perturbation parameter A is in the domain 



'Dxo ■■= {A e 



|A| < Ao}. 



Here Aq is an arbitrary constant, but recall that in case Aq > 1, Assumption 2.7 may not be satisfied for a, b > 1. 
We define the spectral parameter z = E + ir], with E gM. and ?7 > 0. Let -Eq > 3 + Aq and define the domain 



VL--={z = E + i'neC : < ^0 , {<fiN)^ <Nr]< 3N} , 



(2.12) 



with L = L{N), such that L > 12^. Here, we chose Eq bigger than 3 + A, since we know that the spectrum of 
W lies in the set {E e M : < 3} with high probability. Thus spectral perturbation theory implies that the 
spectrum of H is contained injii'eR: !-E|<3 + A}, with high probability. 

Recall the definition of ke, the distance to the endpoints of the support of /x/c, in (2.7). In the following, we 
often abbreviate ke = k. 

2.3.1. Local Laws. 

Theorem 2.9. [Strong Local Law] Let H = \V-\-W, where W satisfies the assumptions in Definition 2.1 and W 
satisfies Assumption 2.2 or Assumption 2.7. Let 



(2.13) 



Then there are constants v > and Ci, depending on the constants Aq, Eq, Aq, 0, Co in (2.4) and the measure /i, 
such that for L > 40^, the events 



Pi \ \m{z) - mfc{z)\ < {^PnY^^ ( min 



Ai/2 



A 



1 



1 

'Nt] 



(2.14) 



and 



n 



both have {^,u)-high probability. 



I lnimfc{z) 1 
Nri Jh) 



(2.15) 



Remark 2.10. If we choose the entries o{V= (vi) to be independent standard Gaussian random variables, we 
have the following result: Under the same assumptions as in Theorem 2.9 (except Assumption 2.2 or 2.7) and 
with similar constants, the events 



fl ||m(z)-m/,(z)|<(<^^)-^«(-^ + ^)| (2.16) 



and 



zeT>L 




both have (^, i^)-high probability. Note, however, that the result only applies to the compact domain Vi of the 
spectral parameter z, for the spectrum of the H = XV + 14^ is not bounded. The proof of the estimates (2.16) 
and (2.17) is similar to the proof of Theorem 2.9 and we refrain from stating it explicitly. 
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For A = 0, we have rufc = nisc, where m^c is the Stieltjes transform of the standard semicircle law. In this 
case stronger estimates have been obtained; see, e.g., [11]. Roughly speaking, in this situation we have the high 
probability bounds 

1 /2 

\m{z) - mse{z)\ < ^ and |G,,(z) - <5,,m(z)| < (^^^^^^) + ^ , (2.18) 

(up to logarithmic corrections), within the range of admitted parameters. 

This suggests that the bound on Gij(z), {i ^ j), in (2.15) is optimal. However, for A 7^ 0, the individual 
diagonal resolvent entries Gn(z) do not concentrate around their mean m(z), due to the fluctuations in the 
random variables {vi). This becomes apparent from Schur's complement formula (see, e.g., (3.16)) and one easily 
establishes that \Gii[z) — to(z)| = 0{\) + o(l), with high probability. 

Comparing the estimate on m — m/c in (2-14) with the corresponding estimate in (2.18), one may sus- 
pect that the leading correction terms in (2.14) stem from fluctuations of the random variables (v,). The next 
theorem asserts that this is indeed true, at least in the bulk of the spectrum: There are random variables 
Co = Co^(^); which only depend on the random variables (wi), but are independent of the random variables (wij), 
such that \m{z) — mfdz) — Co{z)\ < {Nr])~^ with high probability in the bulk of the spectrum; see (2.19). Con- 
cerning the spectral edge, we remark that the estimate in (2.14) is optimal for A -C N~^/^, but it is not known 
whether A^/^TV^^/'' is the optimal rate for A > N^^/^. 

To state our next result, we define the domain 

BL-=T>Lr\{z = E + ir]&<C : ^ke+v > {'Pn)^^N-^/^} . 
Wc have the following result: 

Theorem 2.11. Let H = XV + W , where W satisfies the assumptions in Definition 2.1 and W satisfies 
Assumption 2.2 or 2.7. Then, for any z G Vl, A e T^x^, there exist random variables Co(^) = Co {z), such that, 
with the same constants as in Theorem 2.9, the event 

fl Um{z)-mUz)-Uz)\<{^Mr^^\ (2.19) 

has {^,i')-high probability. The random variables Co{z) have the following property: The event 

Ai/2 A 1 



n |Co(^)|<(rf^«min 



(2.20) 



has {^,v) -high probability. 



Remark 2.12. The estimates in (2.19) and (2.20) need some explanation: Choosing E in the bulk of the 
spectrum, i.e., ke > for some x > 0, we have 

\m{z) - rufciz) - Co{z)\ < M''^^ , \m{z) - mfc{z)\ < (<^jv)'^^ (^-^ + , 

with high probability and the estimate seems to be optimal. In particular, on microscopic scales, rj <^ VV^^/^, the 
local fluctuations stem from the Wigner matrix W, whereas on intermediate scales, rj ~ N~'^/'^, the fluctuations 
due to the Wigner matrix are of the same size as the fluctuations due to the diagonal matrix V . Finally, on 
macroscopic scales, r; ~ 1, the fluctuations are dominated by the matrix V . 



2.3.2. Eigenvector delocalization. Next, let /ii < ... < /i^r denote the eigenvalues oi H = W + W, and let 
Ui, . . . ,Un denote the associated eigenvectors. We use the notation = {ua{i))^i for the vector components. 
All eigenvectors are ^^-normalized. The next theorem asserts that, with high probability, all eigenvectors of 
H = XV + W are completely delocalized: 
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Theorem 2.13. [Eigenvector delocalization] Assume that H = XV+W satisfies the assumptions in Definition 2.1 
and Assumption 2.2 or 2.1. Then there is a constant u > 0, depending on Aq, Eq, Xq, and Co in (2.4) and the 
measure fi, such that for any ^ satisfying (2.11), we have 

max max luaii)] < — , 

l<a<Nl<i<N' "^^1- 

with {£,,u)-high probability. 

Remark 2.14. In case the entries ofV = (vi) are independent Gaussian random variables, the situation is more 

subtle: For any finite Eo, there exists a constant Ceq, independent of A'', and a constant v, depending on Aq, Eq, 
and Co in (2.4), such that for any ^ satisfying (2.11), 

max K(i)|<CB„^^^, (2.21) 

\<1<N y/]\f 



with (^, i/)-high probability. However, ceq oo and — >■ 0, as E'o — >■ oo. 

2.3.3. Density of States. Next, wc state our main results about the local density of states oi H = XV + W: Define 
the normalized eigenvalue counting function of H by 

1 ^ 

p{x)--=j;^T.^^'^-l''^y (2.22) 

a=l 

For El < E2, also define the counting functions 

1 f^^ 
n{Ei,E2) := —\{a : Ei < < -^2}! , nfc{Ei,E2) := / pfc{x)dx, 

-'^ JEi 

where we denote, for consistency of notation, by p/c the density of the free convolution measure /U/c- 

Theorem 2.15. [Local density of states] Let H = XV + W , where W satisfies the assumptions in Definition 2.1 

and, V satisfi,es Assumption 2.2 or 2.7. Let ^ satisfy (2.13). Then there are constants z/ > and C, depending 
on Aq, Eq, Xo, 9 and Co in (2.4) and the measure p, such that for L > 40^, the following holds: For any Ei, 
E2, satisfying —Eq < Ei < E2 < Eq, E2 > Ei + {ipff)^N~^ , and any X e X>Ao; the estimate 

W£„^)-„,.(£..£,)|<to)«(l + J^=1I=_^J , (2^23) 

holds with (^, y)-high probability. 

Moreover, let x > 0. Then, there exists a constant C^, depending only on k, such that for any Ei,E2, 
satisfying Li + k < E\ < E2 < L2 — k, E2 > Ei + {i.pn)^N~^ , and any X G T>\^, the estimate 

\n{Ei,E2) - nf,{Ei,E2)\ < C^imf^ + ^^^^^^^) , (2.24) 

holds with (^, u)-high probability. 

We remark, however, that the estimate in (2.24) deteriorates at the edge: We have Cj^ — )• 00, as x ^ 0. 

2.3.4. Rigidity of eigenvalue spacing. Recall that we denote by /Ui < p.2 < ... < /Ujv, the eigenvalues of 
H = XV + W. The estimates on the density of states in Theorem 2.15 imply the following result on rigidity 

of eigenvalue spacing, which establishes the relation between j/i^ — pj\ and |i — j|. 

Theorem 2.16. Assume that H = XV+W satisfies the assumptions in Definition 2.1 and Assumption 2.2 or 2.7. 
Consider pn < fij in the hulk of the spectrum, in the sense that i > eN and j < (1 — e)N for some constant e > 0. 
Assume that \i — j\ > (i^jv)*^"^; for some constant Co > C, where C is the constant in (2.24). Then, there exist 
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constants C\, C2, depending only on e, uiax.{pfc{x) : fJ^i < x < /ij} and min{p/c(a;) ■ l^i < x < such that the 
following holds: For v > as in Theorem 2.15, the estimate 



N 



N 



holds with v)-high probability, for A e Vx^ . If we assume further that \i — j\ < {ipNY^N^^'^ , for some constant 
c > 0, then there exists a constant K such that 



iMi-Mil 



\^-J\ 



(2.25) 



with {£,,u)-high probability, for all A e ©Ao- 



Remeirk 2.17. The estimate (2.25) can be extended to \i — j\ < (<^jv) '^^N^^'^ in the following sense: There 
exists a constant K such that, for some /x^ e 



iMi - Mil 



NpfM 



1 



<(^iv)"^^, 



with (^, i/)-high probability. The estimate easily follows from the proof of Theorem 2.16. 

2.3.5. Integrated density of states and rigidity of eigenvalues. Define the integrated density of states by 

1 



<E) ■■=^\{a : l^a<E}\. 



Similarly, we set 



nfc{E) :-- 



Pfc{x) dx , 



where p/c denotes the density of the free convolution measure pfc- Then we have the following result: 

Theorem 2.18. Let H = XV + W, where W satisfies the assumptions in Definition 2.1 and V satisfies Assump- 
tion 2.2 or 2.7. Let ^ satisfy (2.13). Then there are constants v > and C, depending on Aq, Eq, Aq, and Co 
in (2.4) and the measure p, such that the event 



fl !^\n{E)-nME)\<{^r,f^(^j: 



A3/2 



A ^ Ay^^E^ 



has {£^,1^) -high probability. 

Our final result concerns the rigidity of the eigenvalue location. We define the 'classical' location of the 
a*^-eigenvalue of iJ, 7c by 



f 

J — ( 



Pfc{x)dx 



1 1 

N 



(2.26) 



where pfc is the density of the free convolution measure pfc- 



Theorem 2.19. Let H = XV + W, where W satisfies the assumptions in Definition 2.1 and V satisfies Assump- 
tion 2.2 or 2.7. Let £ satisfy (2.13). Then there are constants u > and C, depending on Aq, Eq, Aq, and Cq 
in (2.4) and the measure p such that 

\Pa - 7a| < imf^ {n-^^^ [s-V3 + 1 (3 < imfH^ + A^/'7V^/*)] + A^Ar-V3s-2/3 ^ xN-'/^) , (2.27) 

with v)-high probability, for all X G X'ao; where we have abbreviated a := min{a, N — a}. 



9 



Remark 2.20. Let us compare this rigidity result with the corresponding rigidity result for Wigner matrices 
(A = 0): In the bulk of the spectrum, where a ~ A/', we obtain from (2.27), 

K-7.I<M \-^^^) ' 

with i^)-high probability. Thus for A ^ 0, the leading corrections in the rigidity estimate arises from fluctuations 
in the diagonal matrix V , and the eigenvalues do not satisfy as strong a rigidity estimate for their locations as 
in the Wigner case; see, e.g., [18, 19, 15, 16]. However, the eigenvalues satisfy a strong rigidity estimate for their 
relative position or their spacing; see Theorem 2.16 above. For A = 0, the rigidity of eigenvalue spacing is an 
immediate consequence of the rigidity of eigenvalue location. 



3. Weak Deformed Semicircle Law 

In this section, we prove a weaker form of the deformed semicircle law. This weak deformed semicircle law will 
be used to prove the strong deformed law in Theorems 2.9. Moreover, complete delocalization of eigenvectors is 
a direct consequence of the weak law stated in the next theorem. 

Theorem 3.1. [Weak deformed semicircle law] Let H = XV + W satisfy the assumptions in Definition 2.1 and 
Assumption 2.2 or 2.7. Let S, satisfy (2.11). Then there are constants C, v > 0, depending on Aq, Eq, Aq, 0, Cq 
in (2.4) and the measure /i, such that, for 

ao<^<Aolog\ogN, L>12^, 

the event 

fl |max|G,,(^)|<C^^| , (3.1) 

has (^, v)-high probability. 

Denote by , the expectation with respect to the random variable Vi, i G {1, . ■ . ,N}. Then the event 

has {^,v)-high probability. 

Moreover, we have the weak local deformed semicircle law: The event 

n \\m{z)-mUz)\<C-^^]. (3.3) 



zeVL 



has [S^^v) -high probability. 



The rest of the section is devoted to the proof of Theorems 3.1 and 2.13. The proof follows closely the proof 
for Wigner matrices, see, e.g., [12]. We will always assume that W satisfies the assumptions in Definition 2.1 and 
that Ay satisfies Assumption 2.2 or 2.7. 

3.1. Preliminaries. 

3.1.1. Some properties of and rufc. The next lemma collects some useful properties of m/c under Assump- 
tions 2.2 or 2.7. 

Lemma 3.2. There exist Li < L2 such that the free convolution measure fife has support [Li,L2\. The Stieltjes 
transform, m/c, of jifc has the following properties, for all z = E + irj € T>l, A G 2?Ao- 
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i. Let K := min{|£^ — Li\, \E — L2\}, then 

ii. There exist constants C,c> 0, depending on ji, Eq and Xq, such that 

c<\X- z-mfc{z)\<C . (3.5) 

We refer to (3.5) as 'stability bound' and remark that a similar condition has already been used in [36]. The 
proof of this Lemma is given in the Appendix. 

3.1.2. Minors. 

Definition 3.3. LetT c {I, . . . ,N}. Then we define H^'^^ as the {N - |T|) x (A^ - |T|) minor of H obtained by 

removing all columns and rows of H indexed by i €. T. Note that we do not change the names of the indices of 
H when defining H'^'^\ More specifically, we define an operation ttj, i e {1, . . . ,N}, on the probability space by 

{TTi{H))ki := l{k ^ i)l{l ^ i)hki . 

Then, for T c {1, . . . , A''}, we set ttt := IlieT ^''^^ define 

:= {{MHhk.^T. 

The Green functions G^'^\ are defined in an obvious way using Moreover, we use the shorthand notation 

(T) AT 
i i=l 

and abbreviate (i) = {{i}) and, similarly, (Ti) = (T U {i}). Finally, we set 

(T) 



m 



(T) 



N ■ 

Here, we use the normalization N~^, instead {N — |T|)~^, since it is more convenient for our computations. 

3.1.3. Resolvent identities. The next lemma collects the main identities between resolvent matrix elements of H 
and ifW. 

Lemma 3.4. Let H be an N x N matrix. Consider the Green function G{z) = G := {H — z)~^ , z G C+. Then, 
for fee {1, . . . , N}, the following identities hold: 



Schur complement formula: For any i, 



Gii — — ^-^ • (3.6) 

hii - z - hiiGi^hr, 



For i 7^ j. 



For i,j ^ k. 



(y) 

Gij = -GuG^ih,, - J2 hikG^phii) . (3.7) 



G,- = G^ + ^ . (3.E 



11 



Ward identity: For any i, 



^ 1 



(3.9) 



where ri = lmz. 
For a proof we refer to, e.g., [12]. 

3.1.4. Large deviation estimates. We collect here some useful large deviation estimates for random variables with 
slowly decaying moments. 

Lemma 3.5. Let (oj) be centered and independent complex random variables with variance and having subex- 
ponential decay 

F(|ai| > xa) < Ce-^'^', 

for some positive constant C and 9 > 1. Let Ai G C and Bij e C. Then there exist constants ao > 1, Aq > 10 
and Co > 1, depending on 9 and C, such that for ao < ^ < ^o log log A'', and ipN = {logNy° , 



/ N N ^ / ^ \ 

y i=l i=l \i=l J 



l/2^ 



1/2N 



N 



1/2N 



< e 



< e 



< e 



-(logAf)S 



-(log JV)« 



-(logJV)f 



(3.10) 
(3.11) 
(3.12) 



for N sufficiently large. 
For a proof see [21]. 

3.1.5. Schur complement formula. The proof of Theorem 3.1 starts with Schur's formula 

1 



Gii — 



(3.13) 



where, for brevity, Gu = Gu{z), Gjk = Gjk{z), etc. Define Ej to be the partial expectation with respect to the 
i^'^-column/row of W and set 



(i) (i) . 

Zi := (1 - Ei) J2 h^,G^hk^ = Y.^K^^'^^hu - ^4;G«) 
ik 



k,l 
(i) 



(i) 



k kjtl 

here we used hij = Wij + XSijVi. For a family of random variables {Fi,. . . , F^) we introduce the notation 

N 



(3.14) 



(3.15) 



i=l 



Recalling the definition m^'^) = TrG^^^ = J2k^ G^^, wc obtain from Equations (3.13) and (3.14) 

1 



Gi 



Xvi + Wii — z — m(') — Zi 
1 

Xvi- z- mfc - ([v] - yi) ' 



(3.16) 



12 



where 

yi-=Gu-mfc, yi := wu - Zi - {m'^^^ - m) . (3-17) 

Note the difference between Vi and wu: Since we assumed that the (rescaled) entries of W have subexponential 
decay, we have 

K,|<ci^, (3.18) 

with :^)-high probabihty, whereas Vi = 0(1). 

Lemma 3.6. There is a constant C such that, for any z G Vl, A e Vx^ and 1 < i < N, we have 

|m-mW|<-^. (3.19) 

Proof. The claim follows from Cauchy's interlacing property of eigenvalues of H and its minors For a 

detailed proof we refer to [11]. □ 

3.2. A priori estimates on the domain Q{z). Define the z-dependent control quantities 

Ao := max|Gij| , := max , A := |m — m/d . (3.20) 

Note that these quantities also depend on A, but we do not display this dependence, since, as we shall see, 
uniformity in A can always be achieved on the domain Vx^ using the stability bound (3.5). 
Let A e Vxg. For z G Vl, we define an event fl{z) by 

i^{z) := {A„ < (vjv)-'«} n {A < (<^jv)-'«} . (3.21) 

First, we check that we can bound the matrix elements of the Green function of the minor i/^'^ in terms of the 
matrix elements of the Green function of H. 

Lemma 3.7. Let z € T>l, A G Vxg and T c {1, . . . , N}, |T| < 10. Then there are strictly positive constants C, c 
such that the following statements hold with high probability on il{z): 



i. For any i ^T, 



c<|Gi''^|<C. (3.22) 



ii. For any i,j^T, i ^ j, 

cAo<\g\]^\<CAo- (3.23) 

Hi. 

|to_toW|<C'A2. (3.24) 

Moreover, the constants C and c can he chosen uniformly in X € T>Xa . 
Proof. Let z e Vl and A e I?Ao- We will successively use (3.8), i.e., 

- G\f = . (3.25) 

Since we are working on Q{z) we have |Gy| < Ao < (</?jv)~^^, for i ^ j. Next, Equation (3.16) yields 

=\z + TO^'^ - XVi - Wu - Zi\ . 
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By the large deviation estimates (3.10), (3.11), the Ward identity (3.9) and Inequahty (3.19) we have 

1/2 

<C((pAr)«i 



' Im m(') 



my 



(3.26) 



with high probabihty on 0(z). Since A < (i^n) on n{z) and since Nr] > (<y?jv)^^^, we conclude that \Zi\ = o(l), 
with high probability on n{z), for z & Vl and A e I?Ao- Finally, since 



m 



W| 



-21 



on 0(0), by (3.19), we find 



Gi, 



\Xvi - z- mfc\ + 0(1) , 



with high probability on Vl{z). This, together with the stability bound (3.5), proves the lower and upper bound 
on Gii. Note that by (3.5), we can choose the upper and lower bound on Gj, to be uniform in A € 2?Ao> ^ € PL- 
Statements i-lii now follow by iterating (3.25). □ 



Next, we define the control parameter ^'(z), for z € 2?l, by 



(3.27) 



where A = |m — m/c|. Again, we suppress the A-dependence of ^{z) from our notation. We will use \l/ = '^{z) 
to bound various quantities in the following: 

Lemma 3.8. Let z G Dl, A e T>\^. Then there is a constant C such that we have with {£^,u)-high probability 
on Q.{z): 



max \Zi\ < , 

i 

maxID^il < C*. 

i 

The constant C can be chosen uniformly in ^; e ©i,, A e X>Ao- 



(3.28) 
(3.29) 
(3.30) 



Proof. We prove (3.28). Let i ^ j, then by Equation (3.7), the large deviation estimates of Lemma 3.5 and 
Inequality (3.18), 



(y) / , (ii) 

Gij\ < G{\wij\ + kifcGif «;y I) < C(^,v)« + \ ^ E 

k,l \ ' ^'^ 



ilmm^*-') 
Nr] 



with high probability, where we used in the last step the Ward identity (3.9). Since Im^*-'^ — m\ < CA^, by 
Lemma 3.7, we get 



\Gij\<C 



Vn 



(^ 



+ M,{z)]+C^^Ao, 



with high probability. Since lmmfc{z) > Crj, by (3.4), we can absorb the term {ipn)^N ^I"^ into the term ^(2). 
Then taking the maximum over i ^ j, inequality (3.28) follows. The proofs for Zi and are similar. □ 
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3.3. Derivation of the weak self-consistent equation. We now put Equation (3.16) into a form which 
admits an analysis of the average of the diagonal resolvent entries. For n e N, define 

diJ,{v) 



Ruiz) := / j^-^ 
J [\v-z-- 



(3.31) 



m/c(^))" ' 

Note that, by the stability bound (3.5), i?„ is bounded above, uniformly in z and A, for any n. Also note the 
special case Ri = mfc- Recall the definitions [v] = Gu — m/c and |A| = |m — m/c|. 

Lemma 3.9. [Weak self-consistent equation] Let z G Vl. Then there is a constant C such that, for all A € Vx^, 
we have on Q.{z) with {£,,u)-high probability 

A2 



(l-i?2)[v] -RaM \ < C* + C 



log 



(3.32) 



Proof. Since \Xvi — z — nifd is bounded below by (3.5), we can expand Equation (3.16) to second order in ([v] — J^j), 



1 ^ 1 ^ 



i=l 



— ^ Atij - z - mfc 



N 



1 ^ 

-T 



1 



TV ^ 



where, see Equation (3.17), 



^ (\vi - z- mfc) 



yi = Wii ~ Zi~ (m^') - m) . 



N ^ {Xvi - z - mfcY 

{[^]-yif+om-yif). 



(3.33) 



Next, we want to use the 'law of large numbers' to replace the averages in the first two terms on the 
right side of (3.33) by their expectation: It follows from the stability bound (3.5) that the family of func- 
tions gi : T>Xg X 9 (A, z) i— >■ [Xvi — z — mfc{z))~~^ arc jointly Lipschitz continuous with a constant depending 
only on Eq, Aq and /k. Since the (vi) are i.i.d. random variables, McDiarmid's inequality implies, for A S 'Dxg, 
zeVL,n=l,2,3, 



1 ^ 

n12 



1 



" J {Xv-z- 



^ {Xvi - z- mfc) 



mfc 



< Co 



A(yjjv)^ 



(3.34) 



with (^, i^)-high probability. Uniformity in A, z and v can be established by a lattice argument: Choose a lattice 
£ e I^Ao X I^L, with \jC\ < CN^, such that for any (A,^) e 2?Ao x 'Dl there is (A',z') e £, with \z - z'\ < N'"^ 
and |A — A'| < N~'^. Then (3.34) holds for all (A, z) G jC for some sufficiently large Co and some sufficiently small 
v > 0. Using the joint Lipschitz continuity of {gi), we conclude that there is a constant C > Co such that the 
event 

N 



n n 

^=l,2,3(A,z)e•DAoX2'i 



-E 



mfcY' 



Aix(v) 



(Xv — z — mfc)" 



< C 



A((/?Ar)« 



N 



(3.35) 



has (^,z^)-high probability, for some u > 0, depending on Eo, Aq and the distribution /i. 
Hence, 

^Eg,.= / ^ ^^^"^ +R2M+R^[y? + ^j:j. ^3^^ 

^ J ^v-z-mfc N ^ [Xvi - z - mfc)^ 



1 ^ 



^ (Xvi - z- mfc)'- 



:{y!-2[y]yi) + o{{M-yi)^) + o 



/AMI 



with high probability on i^{z), for z G Vl and A G Vxg ■ Recalling the definition of Rn in (3.31) and the functional 
Equation (2.6) for m/c, we obtain 



1 ^ 

(l-R,)[y]=R,[y]^+ 



1 



+ 0(([v]-3^.)^) + o(^), 



1 ^ 



^ {Xvi - z - mfc)' 



:iy!-2{vm 



(3.36) 
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with high probabihty on n{z). Recahing that |[v]| = A, we obtain 

|2[v]3;,|< (^^ + (logAr)niax|3;,|2) , 

(the added factor logiV wih be useful below). Using the estimates in (3.29) and (3.30), Equation (3.33) thus 
becomes 

(l-«*l = «.[vP + 0(i^)+c(^ + ,), 

which holds with high probability on ri(z), z € I?l and A € 2?Ao- Next, observe that, since ImTO/c(z) > Crj, we 
can absorb the third term on the right side of the above equation into the forth term. Finally, we note that we 
can choose the constants uniform in z and A. □ 

To conclude the proof of Theorems 3.1 we reason as follows. Assume, for simplicity, that 1 — R2{z) is bounded 
below (this holds true in the bulk of the spectrum). Recalling that |[v]| = A and the definition of ^{z), we are 
going to show that (3.32) implies 



A < CA^ + O 



( \ 



with high probability on 0,{z). Hence, we obtain the following dichotomy: Either 

for some A/'-independent constant c > 0, with high probability on fl{z), z G T>l, A G T^\o- Using the self- 
consistent equation (3.32), wc establish in the next section, that, for large 77, i.e., ?7 > 2, A + Aq < {lpn)'"^^ , with 
high probability. In other words, ^{z) holds with high probability, for Imz > 2. But then the first inequality 
in (3.37) must hold, for A'' large enough, and we can reject the second inequality in (3.37) for such 77. To extend 
this conclusion to all 77 > {(pp^)^ , we make use of the Lipschitz continuity of the resolvent mapping z 1— >■ G{z), 
which not only allows us to establish that ^{z) holds with high probability for rj small, but also shows that (3.37) 
holds for small 77. This continuity, or bootstrapping, argument is outlined in Section 3.5. This argument applies 
in a straightforward way in the bulk of the spectrum where we have |1 — i?2(-2)| ^ c > 0. For z close to the 
spectral edge, \1 — R2{z)\ can become very small and a slightly modified version of the above dichotomy has to 
be applied (see Lemma 3.12), but the bootstrapping method still applies. 

3.4. Initial estimates for large 77. To get the bootstrapping started, we need estimates on Aq and A, for 
77-1. 

Lemma 3.10. Let 77 > 2. Then for z G Vl, A e Vl, we have 

Ao + A<^, (3.38) 

with {^,i')-high probability. 

Proof. Let X gVl- We &x z G Vl, with rj >2 Then we have the following trivial estimates 

IGlpI < - , Im^l < - , M < - , \Rn\ < f-V , (3.39) 
■> 7] 7] ri \rij 

for any T C {1, . . . , iV}. 

We start with estimating Ag-. From equation (3.7) we obtain using the large deviation estimates in Lemma 3.5, 
that 
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with high probabihty. 

To bound A, we note that 



m < |^i| + |m«-m| + Ki| <C 



with high probability. The self-consistent equation (3.16) can be written as 



1 ^ 



1 



Xvi — z — rufc — ([v] — yi) \vi— z — rufc 



N 
i=l 



Xvi — z — m/c Xv — z — rufc 



The second term on the right side of the above equation is bounded by C with high probabihty, as follows 
from (3.35). To bound the other term, we rewrite it as 



(3.41) 



JV 



-y 

^ ~[ - "^fc - (M - yi)){>^Vi - z- m/c) 

Taking the imaginary part, we see that the denominators of the summands are with high probability larger in 
absolute value than 

for r}>2. Thus, taking the maximum over i, we can bound the right side of (3.41) as 

'(yjv)^' 



3/2 ■ \ y/N J ' 

with high probability. This completes the estimate of A and hence the proof. 



□ 



3.5. Proof of Theorem 3.1. We introduce the control parameters 



a{z) = a := 



l-i?2 



i?3 



p{z) = /3 := 



(7V?7)i/3 ■ 



(3.42) 



Note that for any z S T>l, we have /3 <C {(Pn) Also note that we have chosen (3 to be independent of A. 
Concerning a we have: 

Lemma 3.11. Then there exists a constant K > 1, depending only on Eq, Aq and fx, such that, 



j^VKTrj<a{z)<Ky/]^^, zGVl, A e • 



(3.43) 



In particular, we have lmmfc{z) < Ca{z), for some C > 1 and Rs 1, uniformly in z G Vl and A e T^x^. 

The proof of this lemma is stated in the Appendix; see Lemma A. 7. 

Next, we fix E and vary rj from 2 down to {(Pn)^N~^ . Since y'K + t] is increasing and (3{E + ir]) is decreasing 
in T], we conclude that the equation 



V^^T^ = 2U^KI3(E + i??) 
has a imiquc solution rj = fi{U^ E), for any U > 1. Note that ij{U, E) <^ 1. 



(3.44) 



Lemma 3.12. There exists a constant Uq such that, for any fixed U > Uq, there exists a constant Ci{U), 
depending only on U, such that the following estimates hold for any z G T>l: 



Hz) < Up{z) or K{z) > 
Hz)<C,{U)p{z), 
on 0,{z), with {^,u)-high probability. 



a{z) 
U ' 



ifr,>fj{U,E), 
ifV<mE), 



(3.45) 
(3.46) 
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Proof. Fbi z GVl- Since 



we can write the weak self-consistent equation (3.32) as 

"W = + " (1577) + " + VWAT13^) . (3.47) 



Since ^/|3^A + jS^a < l3y/pK + py/a^ < C{I3'^ + I3a + pA) by Young's inequality, we obtain from (3.47) 

|a[v]-[vp| <o(^^^ +C*{pA + ap + l3^), (3.48) 

with high probability on Q.{z), for some C* > 1. We set Uq := 9(C* + 1). Depending on the size of j3 relative to 
a, we estimate either [v] or [v]^ using the above inequality. We have to consider two cases: 

Case 1: rj > fi{U,E) {a > P, "Bulk estimate") Prom (3.44) we find v'kT?? > 2U'^KP{z) and hence, us- 
ing (3.43) and the definition of C* , 

.a a 

Thus we find from (3.48) with high probability on Q{z) that 

aA < + C*{/3A + a/3 + p^) <2A^ + ^ + 2C*ap . 

Hence, aA < AA"^ + AC'*al3. Thus, we either have aA/2 < AA'^ which implies A > a/8 > a/U (recall that 
U>Uo = 9(C* + 1)), or aA/2 < 4C*a/3 implying A < 8C*/3 < Up. This proves (3.45). 

Case 2: 7] < fi(U,E) (a < /3, "Edge estimate"). From (3.43) and (3.44) we find a < 2U'^K'^p. Thus 
from (3.48), we find 

A^ < 2aA + 2C*(/3A + ap + P^) < C'pA + C p^ , 
for some constant C depending on U. Inequality (3.46) follows. 

□ 

With these lemmas at hand, we are prepared to start the continuity argument: We choose a decreasing 
sequence (%), k = 1, . . . ,fco satisfying fco < CN^, \7]k — %+i| < N~^, r}\ = 2 and r^feg = {(Pn)^N~^. For fixed 
E € [—Eo,Eo] we set Zk = E + ir]k- Recall Lemma 3.12. We fix a U > Uq throughout the remainder of this 
section. 

One easily sees that, for large enough A'', rji > r]{U, E), for any E c [—Eq, Eq]. Therefore Lemma 3.10 implies 
that $l(zi) holds with high probability. This is the starting point of the continuity argument. The next lemma 

extends this result to all k < ko. 

Lemma 3.13. Define the event 

Qk ■■= fi(^fe) n {A{zk) < CC^) {U)p{zk)} , (3.49) 

where 



C^''\U) := 



U ifr,k>fj{U,E), 
Ci{U) ifvk<v{U,E). 



Then, there exists u > 0, such that for any ^, 1 < k < ko 

P(17^) < 3fce-''(l°s^)^ (3.50) 
Note that the estimates in this lemma are uniform in A e I?Ao • 
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Proof. We proceed by induction on k. The case k = 1 has just been proven. Hence, assume that (3.50) holds for 
some k > 2. Then 

nni^,) < n n{zk+i) n ni^,) + n in{zk+i)r) + nni) =:B + a + v{ni) , 

where we set 

A := p([Ofe n {A > ((^iv)-'^}] u [Ofe n {A„ > im)-^^}]) , 
B ■= p(Ofe n ^{zk+i) n {A(^fe+i) > c^^+^\u)p{zk+i)}) . 

We start by estimating A. Using the Lipschitz continuity of the resolvent map z ^ G{z), z G C+, we obtain 

\Gij{zk+i) - Gij{zk)\ < \zk+i - Zk\ sup \G'iAz)\ < N-^ sup — ^ < iV'^ . 

zeT>L zeVL [i^z) 

Thus A{zk+i) < A{zk) + iV^^ < CPizk) < (y^N)-^^ and 

Ao{zk+i) < Ao{zk) + < C^{zk) < (logiV)-2« , 

with high probability on il(zfc), where we used Lemma 3.8. Thus F{A) < 2c~'''^'°s^)^. 

To bound B, suppose first that r]k > v{U,E). Then, using the Lipschitz continuity of the resolvent map we 
find |A(2:fe+i) — A(2fe)| < N~^. Thus we find on flk with high probability 

A{zk+i) < A{zk) + < UP{zk) + iV-6 < '^UI3{zk+i) , 

where we used that /3 is a deterministic decreasing function of rj. 

Suppose next, that rjk > ijk+i > fj{U,E). Then since < all^^, by Equation (3.44), we find, in this case, 

A(zfc+i) < alJ^^. But the dichotomy of equation (3.45) then implies on flk n fl{zk+i) with high probability that 

A{zk+i) < U(3{zk+i). If Vk+i < fi{U,E), the dichotomy immediately yields A{zk+i) < Ul3{zk+i)- This shows 

that B < e-(i°s^)* if r/t > fjiU^E). 

If rjk < viUjE), then also r]k+i < fj{U,E) and hence Equation (3.46) gives A{zk+i) < Ci{U)P{zk+i)- 

Thus, we have proven that, for all k < ko, ^^^k+l) ^ 3e-^('°s JV)« ^ P(O^). This concludes the proof of the 

lemma. □ 

To complete the proof of Theorem 3.1, we need to extend the conclusion of the previous lemma to all z G T>l. 
To accomplish this we use a simple lattice argimient using the regularity of the Green function. 

Corollary 3.14. There exists constants C and u > Q, such that, for ^ satisfying (2.11), 



U ^(^ 



U {A{z)>CP{z)} 



< g-,/(logJV)« _ 



(3.51) 



Proof. We choose a lattice C <Z "Dl with |£| < CN^ such that for any z £ Vl there is a 2:' € £ satisfying 
\z — z'\ < N~^. Using the regularity of the Green function we have for z,z' gVl, 



\G,,iz)-G.,iz')\<r^-'\z-z'\<-. 



N 



(3.52) 



Lemma 3.13 yields 



n |a(.') < f /3(.')} 



z'ec 



> 1 _ g-i'(logJV)f 



(3.53) 
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for some constants C and v. Hence, combining (3.52), (3.53) and N ^ < P{z), we get 



U {A{z)>CI3iz)} 



> 1 _ g-i^(iogAr)« _ 



The first term of (3.51) is estimated in a similar way. □ 
This proves (3.3) of Theorem 3.1. To prove (3.1), we observe that (3.28), (3.42) and (3.51) imply that 



Ao<C 



with high probability on fl{z). Then (3.51) and a similar lattice argument as above yields (3.1). To prove (3.2), 
we note that (3.16) yields 



\E'"^Gu-m\ = 



f dnjv) 1 ^ 1_ 

J Xv — z — mfc{z) N ^ Xvi — z — 



+ 0{[Y] + ma^\yi\), 



mfc{z) N j^^Xvi- z-mfc 
with (^, i')-high probability on ri(z), z G 2?l, A € 2?Ao- From the large deviation estimate in (3.35) we find 



with high probability on fl{z), and we can conclude the proof of (3.2) as above. This finishes the proof of 
Theorem 3.1. 

3.6. Delocalization of eigenvectors. Next, we show that the eigenvectors of H are completely delocalized. 
We denote by Va the normed eigenvector to the eigenvalue of H = XV + W, i.e., 

{XV + W)Va = ^laVa , 

with \\v„\\l = j:,K{i)? hi- 
proof of Lemma 2.13. We follow [12]. For z G 2?l and A G we have 

1 



\Gii{z)\ < 



\Xvi - z-mfc + ([v] - yi)\ 



From the weak deformed semicircle law. Theorems 3.1, we conclude that |[v] — 3^i| = o(l), with high probability. 
Since \vi — z — mfc{z)\ > c > is bounded below uniformly in A G Vx^ and z G Vl, by (3.5), we have 

max|Gii(z)| < C , 

i 

with (^, i/)-high probability, uniformly in ^ G Vl and A G Pao- Set r] = {(fiN)^N~^ , L = 12^. Then, by the 
spectral decomposition of H, 



N 



C > Im Giiina + ijy) = ^ 



(Ma - M/3)^ + 



> 



\Va{i) 



with (^, i/)-high probability. This concludes the proof. 



□ 
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4. Fluctuation Lemma and Strong Deformed Semicircle Law 



In this section, wc prove a fluctuation Lemma (see Lemma 4.1 below) that, when combined with the weak local 
deformed law yields a proof of the strong local deformed law, i.e., Theorem 2.9. Recall that we denote by Ej the 
partial expectation with respect to the i*^-column/row of the matrix W. Set Qi := 1 — Ej. Roughly speaking, 
the fluctuation lemma of this section asserts that, assuming the conclusions of Theorem 3.1, we have 

with high probability, up to logarithmic corrections. For a detailed study of fluctuation averages (for generalized 

Wigner- and band matrices) similar to (4.1) we refer to [14], whose arguments wo follow. The situation for the 
deformed ensembles considered here is in so far different as Qi{Gii) < 0(A), whereas Qi{Gii) <C 1 in the Wigner 
ensemble. Note, however, that Qi{G~^) < {Nr])~^/'^ for the deformed model studied here as well; see below. 

4.1. Fluctuation lemma. Recall the notation A = \m — nifd- We set Qi := 1 — E^, where E^ denotes the 
partial expectation with respect to the i*'^-row/column of the matrix W. 

Lemma 4.1. Suppose ^ satisfies (2.11) and let L > 12^. Let S be an event defined by requiring that the following 
holds on it: There are strictly positive constants C, c, such that 

i. for all z e Vl, A e Vx^, 

k{z) < 7(^) , (4.2) 
where ^ is a deterministic function satisfying ^{z) < (yjv)~^^; 
a. for all z e Vl, A e T>\^, 

Ao < C^{z) < C$(z) , (4.3) 

where 

^{zf := (^^)2, Inim,e(^.)+7(.) 

is a deterministic control parameter; 
Hi. for all z € Vl, A € and any i & {1, . . . , A''}, c < |G,i| < C and 



1 



<c(^^ + ^iz)] <C*(z). (4.4) 



Assume that S holds with {£,,u)-high probability, then there exist constants C,c, independent of X and z, such 
that, for p e N, even and satisfying p < v{\ogN)^~^/'^ , 



p 

< {CpfP {^{z)f^ , (4.5) 

for all z e I?L, A e . 

For the proof of this lemma, we need the following two lemmas: 

Lemma 4.2. Let the event S he defined as in Lemma J^.l. Let ^ satisfy (2.11), and let L > 12^. Then there exists 
a constant C such that, for z e Vl, A e I?Ao; ihe following holds: For any T C {1, . . . , N}, with |T| < (logiV)^"-^, 

max |g(J) -Gu\< G\T\KI , max \Gf \<CK, 

(T) 

on In particular, we have that \Gl^ \ > c, for some c > 0, uniformly in T and z &T>l, A S I'ao- 



E 
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Proof. For / G N we set 



Ti:=ma^{\G^P\ : i,j ^ T ,i j ,\r\ = l} , Ti :=ma^{\ 



^(T') ^ 



: i^r ,\r\ = z}. 



Equation (3.8), i.e., Gij = G^*''' + GikGkj/Gkk, implies that we have on 



Ti<K + CAl «: (ifN) 



In particular, we have on S that jG^f^l > \Gii\ — Fi > \Gii\ — 2Fi > 0, for any k ^ i and z e Vl, A e Pao- 
Assume that there is a constant Co such that {G]- '\ > \Gii\ — 2Fi > Cq for any T' with |T'| < I, i ^ T', and 
z £ T>L, X & T^Xo- Then Equation (3.8) implies 



-2£ 



Ti < GKl < C(Fi)2 < Fi « {^n)-^^ . 



thence 



r;+i < F; + CqF^ , F;+i < F; + CqF^ , 



rm<ri + Co^F2, F;+i<Fi + Co^r2. 



Thus, as long as GqITx < 1/4, we obtain by induction that 

F;+i < 2Fi, f (+1 < Fi + 4CoZ(Fi)2 < 2Fi , 

and|Gi'''^|>Co-\for any i ^ T', |T'| = Z + 1, / < (log A/')^ ^. By induction on I, this proves the desired lemma. 

□ 

Lemma 4.3. Let the event S be defined as in Lemma ^.1. Let ^ satisfy (2.11), with ^ > 3/2, and let 
L > 12^. Assume that S has {^,i')-high probability. Then there is a constant G such that for any p,l £ N, 
with p, I < (log A/')^~^/^, and for any z G Vl, A e T>x^, we have 



E 



< CP . 



(4.6) 



where T c {1, . . . , Af}, with |T| < I, and i ^ T. 



Proof. By Lemma 4.2 wc have |g|P| > c on ^, for any 1 ^ i with |T| < (logA'')^ ^. On the complementary 
event S'^, we use Schur's complement formula (3.6), 



^ (ij) 
— = Awi + Wii - z - ^ WikG 



(iT) 



i4T. 



Then by Cauchy-Schwarz, the trivial bounds \G^^\ < r?"^ < A^, E|/iij|P < A^p and E|Aui|P < Ag, and the 
boundedness of Vl, we find 



E 



G 



(T) 



l(S^) < 



E 



G 



(T) 



2p 



l(S^) 



1/2 



P(S^)V2 < (c- + GAT + GAr3)Pp(S^)V2 < CP , 



where we used that S has (^, i/)-high probability and that p < (log A')^ 

Proof of Lemma 4-1- We illustrate the idea of the proof for the simple case p = 2: 

2 



E 



1 ^ / 1 \ 1 ^ 



i=l 



Gii 



iV2 



i=l 



1 

Gii 



□ 



(4.7) 
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The first term on the right side is bounded by 

N / , \ 2 



iV2 



i=l 



1 



1 



i=l 



1 



(4.8) 



whore wo used that ^ has i^)-high probabihty and that N < C<i'{z), since lmmfc{z) > Crj. To handle the 
second term on the right side of (4.7), we use Equation (3.8) to write 



(4.9) 



for i ^ j. Thus we obtain 

Eg, 



1 

Gii 



where we used that G^*^ is independent of the entries in the i**^-column/row of W, and that, for general random 
variables A = A(W) and B = B{W), E.[{Q,A)B] = E[BEiQiA] = if B is independent of the variables in the 
i*'^-column/row of W. Using Equation (4.9) one more and applying the same reasoning we obtain 



G 



EQi 



1 

Gii 



G 



(i) 
33 



33'^ 33 



1 

Gii 



^33^ jj 



^Q^[-^_]Q3 



c . ./^(')r'. 



GjiGij 



\Q3 



Hence, introducing the indicator function 1(S), 



Using that [G,^p|l(S) < G$, (i 7^ j), \Q,X\ < 2\X\ and \G\-'\ > c on S, this last expression is bounded by G*^. 
It remains to laound the same expression, when introducing the indicator 1(S°). Using Holder's inequality and 
Lemma 4.3 we get 



1 

Gii 



G 



m 



< supE 



Q^ 



Gji Gij 



10, 



^'^33'^jj 



Gij Gji 
.'^33'^ 33 ^ii, 



1(S) 



(4.10) 



E 



Q^ 



1 



1 



33 



< P(S'=)i/2 



Q^ 



1 



1 



G 



33 



< 



0(1) 



(4.11) 



where we used that S is a (^, i/)-high probability event. Combining the estimates (4.8), (4.10) and (4.11), 
Inequality (4.5) follows for p = 2. 

Next, let 4 < p < z^(log N)^~^/'^ be even. Writing p = 2r, we have 



E 



=1^ E ^iio.(<^) n 



il,...,i2r j = l 



j'=r+l 



Gi.,i., 



(4.12) 



For simplicity, we first assume that we can replace the sum over the indices i = {ii,. . . , i2r) by a truncated sum 
where all indices are distinct, i.e., we consider 



i^,...,i2r■ fe=l ^^'"'"^ k'=r+l Vt^^fc'H'/ 



(4.13) 



all distinct 



As in the p = 2 case, we make each factor of Ga in the above expression independent as of many summation 
indices as possible by a expansion procedure that uses the identities 



r^(T)^(T) 



G 



(4.14) 



kk 
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for i,j, k ^T, k and 



/^(T)^(T) 



^(T) ^(Tfc) ^(T)^(Tfe)^(T) ' 
'^ii '-^ii '-^ii '^kk 



(4.15) 



for fc ^ T, fc 7^ i. 

The expansion procedure goes as follows: We start with expanding Fj^ := G^J^ in (4.13). Using formula (4.15), 
where the choice of fc e {ii, . . . , i2r}\{ii} is immaterial, we can add to G~ one upper index k. This results 
in two terms, (Fjji := {G\'^IJ~''' and (i^ijo := -Gi^kGui/Gi^i^GfJ^Gkk- Using formula (4.15) we can further 
expand (F^Ji as (Fjjn + (F^Jio, where (F^Jn = (Gj^^*-j)~\ for k G {ii, . . . ,i2r}\{ii,0 (again the choice of I is 
immaterial), and (Fjj)io is a fraction with two off-diagonal resolvent entries in the numerator and three diagonal 
resolvent entries in the denominator. Similarly, we can split the term (F^Jo = (^ijoo + (-Pii)oi; where we 
applied (4.14) or (4.15) to one resolvent entry of (F^Jo, with an index / 7^ ii, k. Again, there is some arbitrariness 
in the choice of the resolvent entry used for the splitting that can, if desirable, be removed by choosing an ordering 
on the set of all resolvent entries G\^ . We continue the splitting of the terms (Fij)o., hereby generating terms 
indexed by sequences a of zeros and ones. 

The precise procedure is the following. Let Q denotes the set of monomials of resolvent entries of the form 
Gnm, with n 7^ m, T C {ii, . . . ,i2r}\{n,m}, and 1/G„„, T C {ii, . . . ,^2r}\{^^}• Given F G Q, the formulas (4.14) 
and (4.15) define an operation, F i-)- Fi € 0, by adding an upper index, e.g., G^X 1-^ G^Z\ and its complementary 
operation F i-> Fq, e.g., gI™ !->■ G^^-'G^'^/G[.'^\ such that F = (F)o + (i^)i- Composing these operations we 
generate from F = (F)0 e G, elements {F)„ e Q, labeled by binary sequences a. For these operations we use the 
notation a aO and a al. Given F = (F)0 e G, the recursive algorithm is as follows: 

(A) Stopping rules 

(T) 

(1) If all terms in (F)o- are maximally expanded, i.e., each resolvent entry in {F)„ is of the form Gnm with 
n, m ^ T, (Tnm) = {ii, . . . , i2r}', 

(2) else if (F)^ contains 2p off-diagonal resolvent entries in the numerator; 
we stop the expansion. 

(B) Else, we choose an arbitrary resolvent entry Gnm in (-^)(t* If — we use (4.15), with some arbitrary 
k e {ii, . . . , i2r}\{('rn)}, to split (F)cr = {F)„o + {F)ai- H n ^ m, we use (4.14), with some arbitrary 
k e {ii, . . . , j2r}\{(Tnm)}, to split (F)<, = {F)„o + {F)^i. 

Below, we show that the stopping rules ensure that the recursive procedure is terminated after a finite number 
of steps. Choosing F = (Fj^) = G^^, the above procedure yields 



Q^J-P^)■■■Q^P (77^) =EQn(^.J<xQ.J77^) •••0., (77^) (4-16) 



ilii / \'-'«2»2/ \^ipip / ^ J \^ 

where the summation index u runs over a subset of finite binary sequences (the number of terms in the sum is 
estimated below). The summands (Fi^^^ G G arc fractions with off-diagonal entries of G in the numerator (except 
for the maximally expanded leading term (G^^]^)~^) and diagonal resolvent entries in the denominator. All these 
entries are maximally expanded in the summation indices. Each term in the rest term Ri^ , a fraction of resolvent 
entries, contains 2p off-diagonal resolvent entries in the numerator. 

We claim that the total number of terms generated by the above recursive procedure is bounded by {Cp)^^, 
for some p-independent constants G. Indeed, the procedure described above, generates a finite rooted binary 
tree, whose vertices are labeled by the binary sequences a. By the stopping rules (1) and (2), each term on the 
right side of (4.16) corresponds to a leaf node of this tree. Thus to get an upper bound on the number of terms 
in (4.16), it is enough to estimate the depth of the tree. 

To estimate the depth of the tree, we estimate the maximal length of a generated sequences a. We first 
observe that the number of off-diagonal resolvent entries is raised by one under the operation a 1— )■ ctO (except in 
case it is first applied to Fi^, when the number is raised by two). Hence, by stopping rule (2), the leaf nodes are 
labeled by sequences a with at most 2p zeros in it. Also note that the operation a aO increases the number of 
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resolvent entries by at most 4, but the operation a i— >■ al does not change this number. Thus the total number of 
resolvent entries in a term (-Fijo- is bounded by 8p + 1. Hence, a bound on the number of upper indices for any 
vertex on the tree is (8p+ l)p. In other words, a sequence labeling a vertex has at most {8p+ l)p ones in it. Thus 
a sequence labeling a leaf node has a most {8p + l)p ones and 2p zeros, therefore has length at most 8p^ + 3p and 
we conclude that the number of leaf nodes of the tree is bounded by 



for some constant C, independent of p. 

It follows that the right side of (4.16) contains at most [CpY^ terms. In particular, the remainder i?,^ contains 
at most {Cp)'^P terms, each of which contains 2p off-diagonal matrix resolvent entries and less than 3p diagonal 
resolvent entries. By assumptions ii, Lemma 4.2 and Lemma 4.3, the rest term Ri^ thus satisfies 

\Ri\liE) < {CpfP^izfP , E\Rif < {CpY^N^p . 

(Here, we also took into account that each factor of Qi yields a factor 2). The second bound follows from the 
rough bound < r/-^ < N. Thus, for p < v{\ogN)^-^/'^ , 

E|i?,|l(S^) < P(S^)i/2E[|i?,|2]V2 < ^ 

since S has (^, i^)-high probability by assumption. 

Next, we expand the term term C?~ (4.16). We apply the same procedure to each 'leaf node term' 



Qii{Fii)aQi2{GiX) ■ ■ ■ Qip^ijip ■ 

Note that we do not expand the remainder term Ri^ any further nor start a new expansion separately for G~X^ 
(this would yield an expansion with too many terms for our purposes). We also modify the stopping rule (2) 
accordingly: We stop expanding a term in (4.16) whenever it contains 2p off-diagonal resolvent entries. Applying 
the algorithm iA)-{B) to (4.16) we find 



Qh [ 7; )Qi2 [ ) • --QiiT ( -F, 1 — X! QiiiFh)(TiQi2{Fi2)a2 ■ ■ ■ Qi2r ( ) + -Rii + -Rjs ) 

(4.17) 

where the remainder Ri^ satisfies the same bounds as Ri^ . The effect of the modified stopping rule (2) is that 
the sequences cti and together contain in total at most 2p — 1 zeros. 
Expanding the remaining 2r — 2 factors of G^^ in (4.17), we find 



^UQ^A?^) n Q^.'in = E ^ QniF^,U---Qi2AFi2.)a2. +^'^' (4-18) 

where the remainder TZ = Y^g^i Ri^ satisfies 

|7^|1(S) < (Cp)2f ^ ]E|7^|2 < (Cp)4pjv4p . (4.19) 

Thus 

E|7e| < {CpfP^^P . (4.20) 

It therefore suffices to consider only the first term on the right side of (4.18), in which all monomials {Fi^)^^ 
are maximally expanded and the summation runs over 2r binary sequences of finite length. Note that the total 
number of zeros in the array of sequences a = (cri, . . . , a2r) is, by the modified stopping rule (2), at most 2p — 1. 
It follows that the total number of terms in (4.17) is less than {Cp)^^. Indeed, this can be checked in the same 
way as is done above: A term in (4.17) corresponds to a leaf node on a rooted binary tree, whose vertices are 
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labeled by a. The total number of zeros in a indexing a leaf node is bounded by 2p and the number of ones is 
less than (8p + l)p^. It follows that the total number of terms in the expansion of (4.17) is bounded by {Cop)^P 
and we find 



(4.21) 



Recall that, due to our simplification assumption all indices {ii,.. . ,ip) are distinct. As in the case p = 2 we now 
use the presence of the Q's: First, we claim that, for any label a e {1, . . . , 2r}, 



|g,jF,j,ji(s)<(c$)i+°(-»), 



(4.22) 



where O(cto) denotes the number of zeros in the sequence Ua- For 0{ua) = 0, this follows from Hi. If 0{ua) > 1, 
the successive application of the operation cr i-)- aO has generated 0{(Ta) + 1 off-diagonal resolvent entries and at 

most 3 0((To) + 1 diagonal resolvent entries. 

Next, choose (ii, . . . , Z2r) and (di, . . . , (T2r) in (4.21) such that 



E 



The key observation is the following: 



^0. 



(4.23) 



(C) Let a e {1, . . . , 2r}, then there is a label h E {1, . . . , 2r}\{a}, such that the monomial {Fi^)^^, contains an 
oH'-diagonal resolvent entry with ia as a lower index. We use the notation b = 1(a), if b is linked to a in this 
sense. 

Indeed, assuming the contrary, we conclude that all monomials (i^ija^ in (4.23), but (-Fi„)o-„, are independent 
of the random variables indexed by ia- But due to the presence of the Qi^ this term has vanishing expectation. 
Note that this argument relies on the assumptions that all indices {ii, . . . , i2r} are distinct. 

Next, let a e {1, ... , 2r} and denote by la := |l~^({a})|, the number of times the label a is linked to some 
label b in the sense of (C). Then 



|(FiJ<,Jl(S)<Cf$i+'<'. 



(4.24) 



Indeed, for each label c e l^^({a}) we had to use at least once the operation a aO to get the lower index ia- 
Hence, 0{ab), the number of zeros in ah, is at least la. Inequality (4.24) follows from (4.22). 

Finally, noting that la > 1, for each label a € {l,...,2r}, we find that each non- vanishing term of the 
form (4.23) is bounded by (C*)^^. 

Using Lemma 4.2, Lemma 4.3, the bound on the remainder term (4.20) and the assumption that the event E 
has (^, z^)-high probability, we obtain 



»l,...,l2r 

all distinct 



fe=l 



Gi 



It 

n • 

/c'=r+l 



for any p < i^(log A^)^^"^/^, under tlie simplifying assumption that all indcces are distinct in the sum. 

To deal with the general case, we go back to (4.12). Abbreviate i — (ii, . . . ,i2r)- Denote by V2r the set of 
partitions of {1, . . . , 2r}. Let r(i) be the element of V2r defined by the equivalence relation a ~ 6, if and only if 
ia = ib- Then we can write 



E 



nT.Q^[g-) E E i(r = r(i))Eft, 



1=1 



1 



1 



(4.25) 



Fix now i, and denote by F := r(i), the partition induced by the equivalence relation ^. For a label a G {1, . . . , 2r}, 
we denote by [a] the block of a in F. Let S{T) :— {a : |[a]| = 1} C {1, . . . , 2r} denote the set of single labels 
and abbreviate by s |'S'(F)| its cardinality. We denote by is(r) '■— (*a)aeS: the summation indices associated 
with single labels. Notice that if a. is a single label (for some F), then there is exactly one Qi^ on the right side 
of (4.25). However, if a is not a single label (for some F), Qi^ appears more than once on the right side of (4.25). 

Next, we expand the summands on the right side of (4.25), using the recursive procedure {A)-{B), but we 
only expand in the single labels. More precisely, the recursive procedure is now defined by 
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{A) Stopping rules 

(1') If all terms in {F)^^ are maximally expanded in the single labels; a resolvent entry Gnm, is maximally 
expanded in the single labels if ig C (Tnm), n, m ^ T; 

(2) else if {F)cr contains 2p off-diagonal resolvent entries in the numerator; 



we stop the expansion. 



{B') Else, we choose an arbitrary resolvent entry Gnm in {F)^. If n = m, we use (4.14), with some arbitrary 
index k S {is}\{(Tn)}, to split {F)„ — {F)„o + {F)„i. If n ^ m, we use (4.15), with some arbitrary 
k e {is}\{(Tnm)}, to split {F)^ = {FU + {F)^i. 

Applying this procedure, we obtain a similar expansion as in (4.17). The remainder terms can be estimated in 
the same way as before, simply by using the fact that each term in the remainder contains at least 2p off-diagonal 
resolvent entries. Also note that the bounds on the number of terms in the expansion still apply. It therefore 
suffices to bound the summands in (4.18), (now some of the indices may coincide). Recall that s denotes the 
number of single labels in the fixed configuration i. We claim that 



(4.26) 



This follows in a similar way as above, using the following observation: 

(C) Let a € S{T), then there is an label b G {1, . . . , 2r}\{a}, such that the monomial (i^it)<7t contains an 
off-diagonal resolvent entry with ia as a lower index. 

The bound (4.26) now follows in the same way as above, by only considering single labels. 

We now return to the sum in (4.25). We perform the summation by first fixing a partition F e V2r- Then 



1 / 1 / 1 

^i:i(r.r(i,|< - < 



^/N 



2r-s 



(4.27) 



Here we used that any block in the partition T that is not associated to a single label, consists of at least two 
elements. Thus |r| < (2r-hs)/2 = r + s/2. Now, using N''^/^ < C$, we find, combining (4.27), (4.26) and (4.19), 



E 



N ^ 



1 



2r 



1(S) < (CpfP J2 (^*)''' ■ 
reV2r 



Next, we recall that the number of partitions of p elements is bounded by {Cp)^^, thus 



E 



2r 



1(S) < (Cp)5f$2p_ 



On the complementary event E'^, we use as before Lemma 4.2 and 4.3 together with the assumption that the 
event S has (^, z^)-high probability, to obtain a similar bound, for p < {log N)^~^^'^ . We omit the details. □ 



We will use the fluctuation Lemma 4.1 in a slightly generalized setting. Abbreviate 

1 



Xvi — z — mfc{z) ' 



(4.28) 



and also recall that the random variables (5,) are bounded uniformly in A and z as follows form the stability 
bound (3.5). We will use the following corollary of the fluctuation lemma: 



Corollary 4.4. Suppose ^ satisfies (2.11) and let L > 12^. Let E be the event defined in Lemma 4-i and assume 
it has (^, v)-high probability. Then there exists a constant C, independent of X and z, such that, for p e N, even 
and satisfying p < ^{log N)^~^^^, and n— 1,2, 3, 



E 



<{CpfPmz)fP , 



(4.29) 



forzeVL,XeVxo. 
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Proof. In the proof of Lemma 4.1, we used the following two properties of tlie (Qi): 

i. For general random variables A = A{W) and B = B{W), E[{QiA)B] = E[BEiQiA] = 0, if B is independent 
of the variables in the i*'^-column/row of W. 

ii. For a general random variables X, \QiX\ < 2\X\. 

Define Qi := QiQ^. Since the random variables (vi) are independent of the random variables (wij), property i 
holds true with Qi replaced by Qi. (Here E stands for the expectation with respect the (wij) and the (vi) random 
variables, but, since the random variables (gi) are uniformly bounded, one could replace E by the conditional 
expectation with respect the (wij)). Since the family of random variables (gi) is uniformly bounded, property ii 
holds now with \QiX\ < 2\gi\'^\X\ < 2Cg\X\, for some constant Cg, for any random variables X. Thus the proof 
of Lemma 4.1 also applies to left side of (4.29): It suffices to multiply the bounds with C|". 

□ 



4.2. Strong self-consistent equation. With Corollary 4.4 at hand, it is easy to derive a stronger self-consistent 
equation for m — rufc- Recall the notation [Z] = J2f=i ^i- 

Lemma 4.5. Suppose ^ satisfies (2.11) (with ^ > 2). Assume that there exists a deterministic function 7(2;) with 
< {(Pn)^^^, such that for all z S Vl, A € I?Ao» 

A(z) < 7(z) . 

with (^, i')-high probability. Then we have with (^ — 2, v)-high probability 

'lmmfc{z) + 7(z) 



1 ^ / I 



i=l 



with n= 1,2, 3, and 



\[z]\<c{^Nr^ 



Nt] 

Immfcjz) + 7(z) 



(4.30) 



(4.31) 



where the constant C can be chosen uniformly in z and A e I?Ao- 

Moreover, the strong self- consistent equation 



1-R 



i?3 



< O 



A2 



logiV 



+0 i^Nr^ 



Nr] 



+ 



A((pAf)^ 



(4.32) 



holds with (^ — 2, u)-high probability, uniformly in z G Vl, A e Pao- 



Note that in the above lemma we have not changed the value of the parameter u, but replaced the A''-dependent 
parameter ^ by ^ — 2. This is necessary at this point, since in the iteration procedure below we apply this lemma 

log log times. 

Proof. We begin by proving (4.31). From Schur's complement formula we obtain 



Qi 



ii) 



= Qi )<Vi + Wi. 



k,l 



Wii - Qi ^hikG^^}hii 



= Wii - Zi . 

Since \wii\ < ((/?Ar)^7V~^/^, with high probability, we obtain from the large deviation estimate (3.10) 



(4.33) 



1 ^ 
N ^ 



< 



N 
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with i/)-high probability. Hence it suffices to bound the average of the left side of (4.33) to get (4.31). 

Theorem 3.1 and Lemma 3.8 imply that assumptions i and ii of Lemma 4.1 hold with high probability. By 
Lemma 3.7, wc have c < \Gii\ < C, with high-probability. Finally, from the estimate on Zi in (3.29), the 
bound on wa, we conclude that assumption iii of Lemma 4.1, i.e.. Inequality (4.4), holds with high probability. 
Hence the event S, as defined in Lemma 4.1, being the intersection of several (^, i/)-high probability events, has 
(^ — 1/2, z/)-high probability. Thus Lemma 4.1, implies that, for any even p with p < u{logN)^~'^, 



E 



1=1 ^ ^ 



< {CpfP^zfP . 



Applying a high-moment Markov inequality, with p = p{\og N)^~'^, we get (4.31). Note that we have not changed 
the parameter ly here, but have replaced ^ by ^ — 2. 

To derive the self-consistent equation (4.32), we return to (3.36), i.e.. 



1 ^ 

(l-i?2)[v]=i?3[v]2 + -^ 



1 



N {Xvi - z- mfcY 



1 ^ 



fri {^Vi -z- mfc)' 



which holds with (^, z^)-high probability. Recall that, on the event H, 

yi = Wii -Zi- (m(') -m)=wu-Zi + {^^) = Q~ + O ($2) , 

After multiplying with gf, (see (4.28)), and averaging we obtain from Corollary 4.4 and a high-moment Markov 
estimate, 

N 



with (^ — 2, !^)-high probability, that is (4.30). Thus, arguing as in the proof of Lemma 3.9, we obtain 



= O 



A2 



logA^ 



+ OliipN) 



Nri 



O 



Xiip 



N) 



N 



with (^ — 2, i/)-high probability, for any z G Vl, X G 'Dao- 



□ 



4.3. Proof of the strong deformed semicircle law. The proof of Theorem 2.9 is based on an iteration using 
the weak semicircle law, i.e.. Theorem 3.1, and Lemma 4.5. We start with an entirely deterministic lemma: 

Lemma 4.6. Assume that 1 < < ^2- Let < r < 1 and L > 1. Suppose that there is a function j{z) satisfying 



+ 



1-T 



such that A(z) < 7(2), for all z G Vl, A € Vx^ . We also assume that, for z e V^, A e Vx^, 



1 - R 



R3 



"Vl-fvf 



A2 



logiV 



^ a{z) + j{z) 



where a = |(1 — R2)/R3\ was defined in (3.42). Moreover, we assume that A <C 1, ifrj^l. Then 



Hz) < M''^' 



A^ _1_ 

iVi/4 + Nt] 



l-r/2 



(4.34) 



(4.35) 



(4.36) 



for allzGVL, AgPao- 
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//V Nr} J - " ' Nri yml"^ Nr] 



Proof. The proof is based on a dichotomy argument. We set 
Note that ao < 7. Using (4.34) we find 

First consider a < ao: From equation (4.35) we find that 

|[v]P< 
and hence 

(|[v]| - af < + (^.)-^^ + ^)"^ + (^.)^--/^^^^^ 



M ^ — [v] 

^3 



2-T 



«|[v]| < o(l)|[v]r + (¥'iv)^^«^ ( ^ + T^j + + '^IMI ' 



(Note that wc used here (^^v to compensate for various constants, as we shall do below.) Thus taking the square 
root and recalling that a < ao and using the definition of ao, we find 

Ai/2 1 \ 1-^/2 

|[v]|<Cao + (^w)"«^ ' 



iVV4 Nrj 

and the claim follows in this case. 

Next, consider a > ao: Assume first that A < a/2. Then in (4.35) we can absorb the terms [v]-^ and A^/ log N 
into the term a|[v]| and we get 



a\/N Nr] aNrj 

where we used the definitions of 7 and ap. Since wo assumed that a > ao, we get A ^ a/2 if A < a/2. Thus, if 
a > ao, we either have A > a/2 or A ^ a/2. By tlic continuity of A(z) in 77 = Imz, we must have A ^ a, since 
we assume that A(^;) <C 1 = 0{a), for 1. Thus, the claim follows from (4.37). □ 

Proof of Theorem 2.9. We prove (2.14). Let £, = log log A/' and set 

f :=2(loglogiV/log2)+^ 
Note that f < 3^/2 < ^0 log log AT. Let L > 40|! To prove (2.14) it suffices to prove 

n {H.)-rn,.(.)|}<(^.)-«^fmin|l^-^,|5^U^'), (4.38) 



with (^, j^)-high probability. 

The weak semicircle law, i.e.. Theorem 3.1 with ^ replacing ^, gives 

~ / Xi/2 1 \ ~ / \i/2 1 \ ^~2/^ 

A < ( ^TTI + 1^ ) < ' 



A^i/4 ]^r]) -^-^'^^ \^ivi/4 Nr] 



30 



for z e Pi, A G Pao) with i/)-high probability. Thus (4.2) holds with 
Since L > 40^, we also have ^{z) < (t^jv)"^^- Hence, by Lemma 4.5 we have 



1-R 



i?3 



< C 



A2 



log A'' 



Nr] 



for 2; e r>L, A e r>Aoi with — 2,z/)-high probability. Since Imm/c < Ca, by Lemma 3.11, this implies (4.35) 
with ^1 = ^. Also, 7 satisfies (4.34) with ^2 = ^ and r = 2/3. Moreover, since A < 7 < (<^jv)~^^, we have A <^ 1, 
if 7/ ~ 1. Therefore, we can apply Lemma 4.6 with ^1 = ^2 = obtain 



A < (ipN) 



lie 



Ai/2 1 

+ 



1-1/3 



for z e A e r>Ao) with (^ — 2, i/)-high probability. Iterating this process M times, we find that 



1 2 / 1 \' 

~ / Xl/2 1 \ 3UJ 



for z e Vl, X e Vxo, holds with (^ - 2M, i^)-high probability. We choose M = [log log iV/ log 2J - 1, here [-J 
denotes the integer part. Since A^/^AT-V^ + iV-V2 + {Nr])~'^ > cN''^ on Vl, we get 



/ A^ J_ 
V7VV4 + _/v^ 



<c< . 



Thus 



A<(¥'jv)^^«( 



7VV4 ^ NrjJ ' 



(4.39) 



for 2; e "Di,, A e "Dao) with (^ + 2, i/)-higli probability (the factor of 2 comes from the —1 in M). This proves (4.38) 
when 

i^r^r^ A ^ AV2 



In case 



we have 



mm 



[ {'pnY^~^ a a^ 



/2 



> 



and the proof for (4.38) is similar to the above case. Finally, when 



> 



1 / Ai/2 1 
+ 



Nt] - 2 Nt] 



(4.40) 



set 



j{z) := {<pr,y'^ 



A 1 
iVi/4 + iv^^ 
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Then by Corollary 4.5 we have 



a|[v]| < CK' + Ci^^r^ ( ^^^^] + 



with j/)-high probability, where we used that 1mm fc < Ca. Assuming that (4.40) holds, using the definition 
of 7 and (4.40), we have 7(2;) < a{z) and we get, using (4.39), 



((.„)"f-^ + i^). (4.41) 

Hence, combining (4.39) and (4.41) we find, using a simple lattice argument, (4.38). Inequality (2.15) then 
follow from (4.38) combined with (3.28). 

□ 



5. Identifying the leading corrections in the bulk 

In this section, wc identify the leading correction terms to m — mjc stemming from the diagonal random matrix V. 
We define random variables Co{z) = Co{z), which only depends on the random variables {vi), such that, in the 
bulk of the spectrum, the leading correction term in the estimate on \m{z) — mfc{z) — Co| is 0{1/Nr]). Finally, 

we prove Theorem 2.11. 

In this section, wc fix ^ = ^^0+0(1) j^g g^j^^^ choose L > 40^. 

5.1. Preliminaries. Let A = |m — m/d- Recall the definition of i?„ in (3.31). In Lemma 3.11, we showed that 

1 - i?2 ~ y/K + T], i?3 ~ 1, 

for all 2; e 'Dl, a e 'Dao- We will need some more notation. For z gVl, A e ©Aoj n Gf^, set 

.^J_Y^ I [ '^i^i'") 

rn[z)-rn. ^ 2-.^^Xvi-z-mfcY J {\v - z - mf^Y ' 

Recall from (3.35) that |r"„(2)| < {(p^)^-^, with high probability, uniformly in 2: e and A G Pao- Thus, 
combining the above observations, we obtain 

C-^V^+Tj<\l-R2-r2\<C^^+lj, C-^ <\R3+r3\<C, (5.1) 

for a\\ z £ T>L, X G Vx^, with (f, z^)-high probability, for some C > 1. 

5.1.1. Definition of Qq. In order to define Co = Co{z), it is convenient to introduce an event Sq, having (^, i^)-high 
probability, by requiring that (5.1) holds on it. We define Co as the solution to the equation 

(l-ii!2-r2)Co(2) = ri(z) + (i?3 + r3)Co(^)', z&Vl, A € • (5-2) 

First, note that Co is well-defined on Sq. Second, note that Co only depends on (wj), but is independent of the 
random entries {wij) of the Wigncr matrix W. Third, from the discussion in the preceding subsection, we infer: 

Lemma 5.1. There is a constant c, such that, for z & Vl, A € I'ao) '"'e have on Sq, 

r A^/^ A 1 

ICol < ((/^at)''^ min 



We omit the proof and just remark that it is suffices to consider the cases K + rj |1 — i?2| ^ ^3 and 
r3< |l-it!2p. 
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Recall the (strong) self-consistent equation for m{z) — mfc{z) in (4.32). The definition of is natural in the 
sense that it embodies the leading correction to m — m/c stemming from the random matrix V: Subtracting the 
defining equation for <^o from the self-consistent equation (4.32), we obtain, after some manipulations. 



(1 - i?2 - r2)(m - m/e - Co) = (^3 + rz){m - mj^f - {R3 + r^Ko + O(A^) + O (^at) 



Imm/c + A 

mi 

on some high probability event S. Theorem 2.11, now follows easily from analyzing the stability of this equation 

in the variable ({z) := m{z) — m/c(z) — Co{z)- 

5.2. Proof of Theorem 2.11. Next, we carry out the details of the proof of Theorem 2.11, which were outlined 

in the previous subsection. 



Proof of Theorem 2.11. Let 

7(^) := {^nY^^ 



mm 



A 



Ai/2 



1 

'Nt] 



Z&T>L. 



(5.3) 



Recall the event Sq, as defined in (5.1). Choosing ci sufficiently large in (5.3), we can achieve that |Co| < l{z) 
on Sq. Next, it follows from Theorem 2.9, Lemma 4.5 and Lemma 3.8, that there is an event Si, having (^, z/)-high 
probability, such that the following holds on it: There is a constant cq such that |A(2;)| < 7(2;), 



max 13^, 



(5.4) 



and, recalling (4.28), 



^|:.."y.(.)|<(.,.r^( '°"''^'^^^+^'" ) 



(5.5) 



for n= 1,2,3, where = wu — Zi — (m*^*-* — m); see (3.17). Since both events Sg and Si have high probability, the 
event S := Sq fl Si has (^,z^)-high probability, with a slightly smaller v > Q. Set (^{z) := m{z) — nifdz) — Co{z). 
Subtracting the defining equation of Co, from equation (3.33), we obtain, using the bounds (5.4) and (5.5), the 
equation 



(1 - i?2 - r2)C = {R3 + r3)(m - m/e)' - (i?3 + r3)Co + ^^(A-') + O {^n) 



^co^ 



Nri 



(5.6) 



on S, for z G Vl, A e ©Ao- Let C2 > max{co, ci} and set 

AV2 



ao{z) := {tpNT^^ I min 



+ 



Thus, on S, we have A <C ao and |Co| ^ cto, for all z G Vl, A € I>Ao- 
Recall that we defined the domain 

BL=VLn{z = E + ir,€C : > M^^^N-^/^} . 

On the domain Bl, it is easy to check that 



mm 



A 



Ai/2 



A 



First, consider the case 



1 1 

> 
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In this case, we can easily see that 



a 



l-R2-r2 



> ao , 



on ^0 n Si. Then we obtain from (5.6), 

m - rufc + Co 



|m - m/c - Col < 



m - m/c - C,q\+U ho — 

a a Nrj 



< o(l)|m - TO/c - Col + O h G- 



a 



on So n Si. Recalling that a + r/ ~ Im m/c on So, we obtain, for some c > C2, 

|m-m/c-Co|<(v'^r« ((^T^A^ + ]^) ' 

on So n Si, for all A e T>\^. Since the condition \/k^\^ > {ipN)^^N~^/^ implies that 

A3 1 11 

< ^TT r < 



we conclude that 



on So n Si , for A € ■ If 



we may use the bound 



{k + 7V3/2 n{k + 7]) - Nr]' 



\m - rufc - Col < 



Nr] 



A 1 1 

< 



|m- m/c -Col < A+ICol < 27 < 



Nt] 



on S, for A G Vx^ . This finishes the proof. 



□ 



6. Density of states 

In this section, we prove Theorems 2.15, 2.16, 2.18, and 2.19. Recall that wc denote by the eigenvalues of 
H = XV + W. In Section 2.3 we introduced the following notations: For Im z > 0, 

N 



and, for Ei < E2, 



n(^i,^2) = ^\{a : ^1 < Ma < E2}\ , n{E) = ^\{a : < E}\ 



Similarly, we have denoted 

fE2 



/Of (x) dx f 
-^^^ — , nfc{Ei,E2) = / p/c(a;)da;, n/c(-E) = / pf, 
X - Z J El J -00 



{x) dx , 



where p/c is the density of the free convolution measure ^/c- Throughout this section we fix ^ = log log A'' 

and choose L > 40^. 
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6.1. Local density of states. Recall that ke ■= min{|ii; — Li\ , i = 1,2}. In the following, we set rj := N~^. 
The first part of Theorem 2.15, Inequality (2.23), is an immediate consequence of the next two lemmas. Their 
proofs follow closely the proof of Lemma 8.1 and Lemma 8.2 in [12]. 

Lemma 6.1. Let rj := N^^. For any Ei < E2 in [—Eq,Eo], we define f{x) = fE1.E2.nix) to be an indicator 
function of the interval [Ei,E2], smoothed out on a scale rj, i.e., f{x) = 1, for x € [Ei^E^], f{x) = 0, for x in 
[El -7], £2 + 7]]", \f'ix)\ < Crj-^ and \ f"{x)\ < Crj-'^ . Assume that the event 



Ai/2 



A 



+ 



Nr) 



(6.1) 



holds with (^, u)-high probability for L := Cq^, for some constant Co > 0. Abbreviate 

K := mm{KEi,KE2} , £ ■= max{£2 - Ei, {ipn)^N~''^} . 

Then, we have 



j f{x){p- Pfc){x)da 



(6.2) 



with {^,L')-high probability. 
Proof. For convenience denote 



p'^ := p- Pfc, m/^ := m - m/c 



We apply the Helffer-Sjostrand functional calculus. We set yo ■= {(Pn)^N ^ and choose a smooth cut-off function 
X such that: 



X{y) = l, on [-£,£]; xiv) = , on [-2£,2£Y; \x'{y)\<j. 

Starting from the Helffer-Sjostrand formula, 

^yf"{x)x{y) + i{f{x) + iyf'{x))x'{y) 



X — X — iy 



dx dy , 



(6.3) 



(6.4) 



we obtain 



J /(A)p^(A)dA <cj dxj^^dy{\f{x)\ + \y\\f'{x)\)\x'{y)\\m^{x + iy)\ 



+ C 



J ^x j\yf" 



{x)x{y)y^^m^{x + iy) 



+ c 



dy/"(a;)x(y)ylmm^(x + iy) 



(6.5) 



Using that x' is supported on [£, 2£], we can bound the first term on the right side of the above inequality by 



■2£ 



dx dy{\f{x)\+y\f'{x)\) 



A 



y/K^^/N SN) 



£X 



1 

+ T7 



(6.6) 



with (^, z^)-high probability. In order to bound the two remaining terms in (6.5), we first bound the imaginary 
part of m^{x + iy). For y > yo, we can use (6.1). So assume that < y < yo- Using the spectral decomposition 
of AU + W, it is easy to see that the function y ylmm{x + iy) is monotone increasing. Thus 



ylmm{x + iy) < yoImm(a; + it/o) < yoImm/c(a; + iyo) + {VN)'^^yo 



Ai/2 



,iVi/4 Nyo 

Recalling that, by Lemma 3.2, we have Imrn/c(a; + iy) < C^/kx + y, we get 

■ Ai/2 1 



{y < yo) . (6.7) 



ylmm(a; + iy) < yoC^/n^ + y + {ipN)^^yo 



iVV4 ^ Nyo 



iy<yo), (6.8) 
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with i/)-high probability. Using that y < yo = {'fiN)^N ^, we can now easily bound 



\ylmm {x + iy)\ < 



N 



iy < yo) : 



(6.9) 



with i/)-high probability. Since by assumption we have 77 < we can bound the second term on the right 
side of (6.5) by 



N 



j <ix\f"{x)\ £dyx{y) 



< 



N 



(6.10) 



with (^, i/)-high probability, where we used that the support of /" has measure 0(77). To bound the third term 
on the right side of (6.5), we integrate by parts, first in x then in y to find the bound 



C 



J dxf'{x)r]Ilem'^{x + ir]) 



C 



da;/ dyf{x)x{y)yRem {x + iy) 



+ C 



J j dyf'{x)x{y)'Rem^{x + iy) 



The second term in (6.11) can be bounded similarly to the first term of (6.5) and we obtain 



dx / dyf\x)x{y)yRem'^{x + iy) 



< ifN) 



with (^, i/)-high probability. To bound the first and the third term in (6.11), we write, for y < yo, 

rvo 

|m^(a; + iy)! < |m^(a; + iyo)| + / du{\dum{x + m)\ + \dumfc{x + m)\) . 



(6.11) 



(6.12) 



(6.13) 



The first term on the right side of (6.13) can be estimated using (6.1). For the others we observe that the Ward 
identity (3.9) implies, for u < yo, 

\dum{x + m)\ = |-^ TrG^(a; + m)\ < X! + — ~ I™'Ti(a; + iu) < -^yolmm{x + iyo) • 



Similarly, we obtain 

.„ , . f pfc{t)dt 1 ^ , . , 1 ^ , . , 

\Oumfc{x + m)\ < / T— — — — 2 = - lmTO/c(a; + m) < — 2/olmTO/c(x + lyo) • 

J t X Itl XL XL 

Prom (6.13) we hence obtain 

\m'^{x + iy)\ < {<pNfi (1 + r du^) < {^Nf^^ , {y < yo) , 
with (^, i/)-high probability. Thus we can bound the first term on the right side of (6.11) by 



j dxf'(x)r}'RBm^{x + ivi) 



< 



(VAf) 



N 



(6.14) 



(6.15) 



with high probability. To bound the third term on the right side of (6.11), we split the integration in the y 
variable into the pieces [?7, yo) and [yo,oc). Using (6.14) we can bound the first piece by 



//•yo 
dx\f{x)\j dy|m^(x + iy)| < 



(V'iv) 



N 
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with high probabihty. For the second integration piece we find 



A 



/d.|/'(x)|£ d„|™-(x + i„)|<to)«/d.|/'(x)|£ d,(^^ , ^.^ 



'yo 

dy , 

'1/0 Vi^ + y, 



< (^iv)^^ ( ^ 



1 



with high probability. Adding all the contributions together, we have proven that 

/(A)p^(A)dA 



with i^)-high probability. 

As a simple corollary, wc obtain: 
Lemma 6.2. Under the assumptions of Lemma 6.1, we have for 



□ 



(n(i;2) - n{Er)) - - n},{E,))\ < {^^f^ ( - + 



1 



(6.16) 



wit/i {£,,u)-high probability. 
Proof. Observe that 

\n{x + rj) — n{x — 77)! < Cr]lm.m{x + irj) < 
with (^, !/)-high probability, where we used (6.7). Hence, 



N 



n(£i) - n{E2) - J /(A)p(A)dA < C {n{Ei + 77) - n{Ei - ??)) < 



N 



(6.17) 



with (^, !^)-high probability. Since p/c is a bounded function, we find 

nfciEi)-nf,iE2)- J /(A)/9/c(A)dA 

Combination with the claims of Lemma 6.1 yields the statements. 

The first statement of Theorem 2.15, i.e., (2.23), now follows easily from the two preceding lemmas. 



□ 



6.2. Bulk fluctuations. The aim of this section is to prove the second part of Theorem 2.15, i.e.. Inequal- 
ity (2.24). Recall the definition of the random variables Co in (5.2). Since, we will restrict the discussion to 
the bulk of the spectrum, where k > 0, we may use slightly modified random variables, Co (-2) = Co{z), which 
approximates Co in the bulk, that are easier to handle in computations. 

6.2.1. Definition o/Co- We define a random variable Co(-z) = Coi^) by 



(6.18) 



for z € Vl, a e T>Xa, where Rn has been defined in (3.31). Recall that, 1 — R2{z) ~ y/n + r) and Rs{z) ~ 1. 
Hence, by the large deviation estimate (3.35), 



ICo|< 



(^iv)"^A 



(6.19) 



with (f , i/)-high probability, for some c, uniformly in z and A. Also note that Co approximates the random variables 
Co in the bulk: Since 1 — i?2 ^ 1 in the bulk, it is straightforward to show that |Co('2) — Co('2)| = 0{N~^), with 
high probability for such z. 
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Lemma 6.3. Under the assumptions of Theorem 2.15, there is c > 0, such that the event 



Pi < \m{z) - rufciz) - C,o{z)\ < {(Pn)'^ ( mii 



zeT>L 

has v)-high probability. 



A 



A2 



+ 



Nri 



(6.20) 



We omit the proof of this lemma since it is similar to the proof of Theorem 2.11. Note however, that the 
estimate in (6.20), deteriorates at the spectral edge, i.e., as k — ^ 0, and we have to restrict the discussion below 
mostly to the bulk of the spectrum. _ 

Next, wc show that the random variable Qo{z), z = E + iri, is a slowly varying function of E (for fixed r]), in 
the bulk of the spectrum. 

Lemma 6.4. Under the assumptions of Theorem 2.15, the event 



n 



9Co(^ + i??) 



dE 



< 



(k + 77)3/2 ViV j ' 



has {£^,y) -high-probability. 

Proof. Recalling the definition of Co in (6.18), we compute, for z G Vl, A e T^x^, 



dE 



d 



1 ^ 



1 



dEl-R2{z)J \Nf^^^"'Xvi-z-mfc{z) 
1 + m'f^jE + i7?) / 1 



+ 



1 - R2{Z) 



N 



^ (An 



{Xvi - z- mfc{z)y 



where we abbreviate m'j^{E + ir/) = '^"'■^ g^^"''' and Qy. := 1 — E^. , where E^. denotes the partial expectation 
with respect the random variable Vi. Differentiating the functional equation (2.6) for m/c, we get 



dlJ,{v) 



mfc{E + iri))' 



(l + m'^,(£ + ir,)) , 



hence. 



K 



for some constant K > 1, where we used Lemma 3.2. Similarly, 

d 1 _ 2(1 + m'f^{E + IT])) f dii{v) 



I {Xv 



dEl- R2{E + i'n) [l- R2{E + \'n)Y J [Xv- z-mfc{E + iri)f 
_ 2(1 + m'f^{E + iri))Ri{E + irj) 



{l-R2{E + \n)f 



and hence, by Lemma 3.2, 



d 



dEl-R2{E + ir]) 



< 



C 



{k + ?7)3/2 



for some constant C. Next, we bound the terms involving the Qy^ by the large deviation estimates (3.35). Finally, 
uniformity in A and z follows from a lattice argument using the stability bound (3.5). □ 
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Next, let /(x) = fEi,E2,r}{x) be an indicator function of the interval [.Ei, £'2], smoothed out on scale rj = N~^. 
Let x{y) be a smooth cut-off function as defined in (6.3). We set m'^ := m — rufc- Appealing to the discussion 
in Section 6.1, we define 

Xo{E,,E2) ■.= ^ [ Axdy {iyf"{x)x{y) + i(/(x) + iyf'{x))x'{y)) Co{x + iy) • (6.21) 

Lemma 6.5. There is a constant C, such that for Ei < E2, with Ei,E2 & [—Eq, Eq], and for A e 'Dxo, we have 

with {£,,u)-high probability, where £ = max{£;2 — Ei, (ipN)^ and k = mm.{\Ei — Li\ : i = 1,2}. 
Choosing the energies Ei, E2, such that mm{KEi, KE2} ^ ^> for some x > 0, we obtain 



6.22 



\Xo{E,,E2)\<CUmf^^, 

with (^, i/)-high probability, for some constant Cj^, depending on x. 
Proof. Starting from the definition of Coj we find 



(6.23) 



\Xo{EuE2)\<C 



da: / dy f{x)x'{y)Co{x + iy) 



/ dyyf'{x)x'{y)Co{x + iy) 



(6.24) 



+ C 



j dx j dy/"(a;)x(y)yImCo(a; + iy) 



To bound the first term on the right side of (6.24) we integrate by part in the variable y to find, with (^, z/)-high 

probability. 



J dxfia 



2£ 



) / dyx'{y)Co{x + iy) 



J dxf{x) dy x{y)dyCo{x + iy) 
j dxf{x) dy x{y)dxC,o{x + w) 



I 



dxf{x) 



{K,+£fl^^/N 



£X 



1 



{K + £)3/^y/N 



(6.25) 



where wc used in the second line that Co{z) in an analytic function in the upper half plane, and in the third 
line we used Lemma 6.4. In the fourth line we used that Kx > k = mhi{KEi, KEi}- Finally, we used that / is 
supported on [Ei — r],E2 + rj], rj = N~^. 

To bound the second term on the right side of (6.24), we integrate by part in the variable x and find, similarly 

to the computation above. 



E2+V 



Ei-ri 



dx / dyyf'{x)x'{yKo{x + iy) 





pE2+rj p 








JEi-T] Jo 



< 



dx / dyyf{x)x'{y)da;Co{x + iy) 
Jo 

{^Nf^ / dx / dyyf{x)\x'{y)\ 

JEi-ri JS 



\ 1 fE2+ri pzc 



r2£ 



£^X 



1 



{K + £)V^y/N 



(6.26) 
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with i/)-high probability. 

Finally, the third term in (6.24) can be bounded by integrating by parts in x to obtain 



f"{x)x{y)y Im Co {x + iy) 



J j dyf'{x)x{y)yir^da;Co{x + iy) 
Jo (k 



£^X 1 



(« + £:)3/2^' 

with i/)-high probability. Adding up the estimates (6.25), (6.26) and (6.27) yields the claim. 



(6.27) 
□ 



6.2.2. Local eigenvalue density in the bulk. In this subsection, we show that we can control the diflFerence 
n(i?2) — n{Ei) in terms of nfc{E2) — UfdEi) in the bulk of the spectrum up to an optimal error: Fix some h > 0. 
We consider energies Ei < E2, such that mm{KEi, KE2} > >c, Li < Ei < E2 < L2 and E2 — E\> {(Pn)^^N~^. 
We denote with constants that only depend on k (with — >■ 00, as >^ — )■ 0). 

As above, let f{x) = fEi,E2,ri{x) be an indicator function of the interval [Ei,E2\, smoothed out on scale 
T) = N~^. Let x{y) be a smooth cut-off function as defined in (6.3) and let := m — m/c. Define 



3ii{Ei,E2) 



2tt 



and recall the definition of Xq in (6.21), 



3Co{Ei,E2) 



1 

2^ 



iwf"{x)x{y) + i{f{x) + Wf'{x))x!{y)) m^(a; + iy) , 



i^y^"{x)x{y) + V{x) + iyfix))x'iy)) Co{x + iy) . 



(6.28) 



(6.29) 



Here we implicitly assume that the functions / and x in both definitions agree. 
Following the discussion in Section 6.1, one easily sees that 



{n{Ei,E2)-nfciEuE2))-Xi{Ei,E2)\ < 



N 



(6.30) 



with (^, i/)-high probability. Recalling the estimate on 3Co in (6.22), we see that it suffices to bound Xi — Xq in 
order to control the density of states. 



Lemma 6.6. Let Li < Ei < E2 < L2, with min{KEi, KE2} ^ ^ E2 — Ei > {ipN)^^N ^. Then 

(£1, £2) - 3eo(^l,^2)| < C^iifNf^^ , 



(6.31) 



with {S,,v)-high probability. The constant C can be chosen independent of Ei,E2 and A € T>\g. 

Proof. We set yn := {(Pn)'"N~^ and abbreviate C,{z) = m{z) — mfc(z) — C,o{z). Using the definition of Co) we find 

l(Xi-Xo)(£;i,i?2)| 



<C j Ax 



{\f{x)\ + \ynx)\)\x:{y)mx + iy)\ + C 



<iyf"{x)x{y)y^T^rn (x + iy) 



C 



d2//"(a;)x(y)y ImCo(a; + iy) 



da; / dyf"{x)x{y)yli^C,ix + \y) 



3.32) 



Using (6.20) we can bound the first term on the right side of (6.32) as 



jdxf\\f{x)\ + y\nx)\)\x!{y)\ 



\C{x + iy)\ 



<C.(^^)^«-, 



- Idxj\y{\f{x)\+y\f'{x)\)^ 



(6.33) 
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with i/)-high probability. 

The second term on the right side of (6.32) is, by (6.19), bounded by 



/rVo ^ ^ rVo ^ 

dx dyf"{x)x{y)ylTXiCo{x + iy) < {ipNf^- dyx{y)y-^===-j= 
Jo vJo vi^ + y VN 

1 rvo \ fTi 



V Jo vN 



< 



N 



(6.34) 



with (^, !/)-high probabihty. 

To control the third term, we note that both functions y ^ y Imm(a; -\-iy) ,y Imm/c(x + iy), are monotone 
increasing. Thus we get from (6.1), 



y 



Imm(a; + iy) < yolm.m{x + lyo) < (Viv)'^^yo ( \/« + 2/o + -^^^ + 



Nyo 



N 



[y < yo) , 



and 



Since yo = {(Pn)^N ^, this yields 



/Imm/c(a; + iy) < yo Imm/c(x + iyo) < + yo , {y < yo) ■ 



y|ImC(a; + iy)|<^^^, {y<yo), 



(6.35) 



< 



(VJv) 



with (^, i/)-high probability. The third term on the right side of (6.32) is thus bounded as 

^^E^jdx\nx)\ f\yx{y) 

with (^, i/)-high probability. 

^ ' - - - - integrate first by parts in the variable x and then in y, to find the 



/rvo 
dx J dyf"{x)x{y)ylmm^{x + iy) 



N 



(6.36) 



To bound the fourth term in (6.32), we 
bound 

p2e 



+ 



'ixf'{x)x{yo)yo R-e C{x + iyo) 



( dx I dyf'{x)dy{x{y)y)'ReC{x + iy) 
•I Jyo 

Using the a priori high probability bounds 

/ xi/2 ^ \ „ 
\m^{x + iyo)\ < {<fiNf^ + ^ j < i^Nf^ , \Co{x + iyo)\ < M 

we bound the second term on the right side of (6.37) as 

J dxf'{x)yoC{x + ir]) 

with (^, i/)-high probability. Hence it remains 



(6.37) 



A 1 



V'^ + yo Vn 



to bound the first term on the right side of (6.37), 
dx dyf{x)dy{x{y)y)'ReC{x + iy) < dx dyf'{x)x'{y)y'ReC{x + iy) 

J Jyo J Jya 

j dx j dy/'(a;)x(y)ReC(a; + iy) 
J Jyo 



+ 



For the first term on the right side, we use (6.20) to find 

dx I dyf'{x)x'{y)y'ReC{x + iy) 

J Vn 



(6.38) 
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with i/)-high probability. Using once more (6.20), we bound the second term on the right side of (6.38) as 

dx / dyf'{x)x{y) Re ax + iy) < C^{^Nf^ / dy— < {^n^^^ , (6.39) 

with (^, i/)-high probabihty. Adding up the different contributions, we find (6.31). □ 

To conclude this subsection, we prove (2.24) of Theorem 2.15: 
Proof of (2.24). Let Ei < E2. Then we have from (6.30) 

\n{Ei,E2)-nfc{Ei,E2)\ < \Xi{Eu E2)\ + CiifNT^ ^ 

< \Xo{EuE2)\ + \Xi{Ei,E2)-Xo{Ei,E2)\+C{ipNr^j^, 
with (^, i/)-high probability. Using Lemma 6.5 and Lemma 6.6, we therefore get 

\n{EuE2) - nfc{EuE2)\ < C^im)"^ f 4 + 



with (^, i/)-high probability. Inequality (2.24) follows by choosing E2 — Ei > {(pn)^^N^^. □ 

6.3. Eigenvalue spacing in the bulk. In this subsection, we prove Theorem 2.16. 
Proof of Theorem 2.16. Let A e T>\^. Starting from the identity 

i-j 



N 

we obtain from (2.24) that 



= Xl{lli) - Xl{llj) , 



with (^, iv)-high probability, for some C large enough. Then, using Ufd^ii) —nfcifJ-j) = (w — Mj)'^/c(K)' some 

where wc used that n'f^{ii^) = Pfdf-i'i) > and 1/C" < p/c < C" in the bulk for some constant C" > 1, depending 
on A and jj. Since — /Ltj| = 0(1), we have 



N 

which shows that the second term in the right side of (6.40) can be absorbed into the left side. Similarly, the 
last term on the right side can be absorbed into the first term in the right side, as we can see from the condition 
> {^Nf^- Thus, 

Ci^<|Mi-/x,|<C2^, (6.41) 

with (^, j/)-high probability for some constants Ci, C2 > 0. This proves the first part of the theorem. 
If \i - j\ < {ifiNf^N^''^, we find from (6.41) that l^i - Mj I < C2{'PnY^N-^/'^ . In this case, 
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with high probabihty. Furthermore, since — < — Hj \ and pfc is analytic (see Remark A. 3), we get 

\pM) - PfM\ < M""^^ , (6.42) 

hence 



N pfcip-dpfM 



with (^, i/)-high probability for some constant K. Thus, we obtain that 



with (^, i/)-high probability, proving the second part the theorem 



< (^jv)^«iV-i 



□ 



6.4. Integrated density of states. The goal of this subsection is to prove Theorem 2.18. The proof follows 

closely [12]. 

6.4.1. Estimate on \\H\\. As a first step, we need an estimate on the operator norm oi H = XV + W. We have 
the following result: 

Lemma 6.7. There is a constant C, such that for all A e ©Aoj we have 

||i?||<max{|Li|,L2} + (^ivf« (^^ + A^) ' 

with {^,L')-high probability. 

Proof. We will only consider the largest eigenvalue /ijv- A bound on the lowest eigenvalue /xi is obtained in a 
similar way. Prom the strong local law (2.14), we get 

Hz) < {^nT^ + ^) ^ z&Vl, a e Pao , 



^Ari/4 TV?? 

with (^ + 2, i/)-high probability. Then we can apply Lemma 4.6, with 

Ai/2 1 



:= {ipN) 



to get, for some sufiiciently large constant ci. 



l-R 



< C 



A2 



r(,n Y^i ( ^ ^ lvamfc{z)+-i{z) 



with (^, i/)-high probability, for any z G Vl and A G Dao- Now, \i E > L2 and k> r],we have lmmfc{z) ^ v/V^ 
and 

I-R2 r- 

a := — V«, 

R3 

by the Lemmas 3.2 and 3.11. Thus, we obtain, upon using Young's inequality. 



R:^ 



log AT 



1 1 
+ 



Vat (Nrj)^ 



(6.43) 



with (^, i/)-high probability, for some ci sufficiently large. 

Given ci, it is straightforward to check that there is a constant C2 > 2ci, such that, for any E satisfying 



L-2 + {vnT^^ 



A 



1 



< £■ < iJo , 



(6.44) 
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we have 



mm{N-^'^K^'\ 7V-i/2a-i«V2^ > (<^^)ci«+2 1 _ (6.45) 



We assume now that E satisfies (6.44) and set 



Note that z = E + ir] G Vl. Prom (6.45), we have k> r]. Similarly, we have 



Imm,.(E + i,)~-^«^; ^«^. (6.46) 



Furthermore, since a > ^JUjK^ for some if > 1, we must have 

A(.)<c(^;v)-^+i(^ + i-)<a, (6.47) 

with (^, i/)-high probability. Here, we used that = (<pjv)'''*+^ and y/R > {ipNY'^^/'^X^/'^N-'^/'^. Since a > A, 

we get from (6.43), 

A < Ci^pNr^ + + , (6.48) 

\a\/N a{Nr]y aN^J 

with (^, i/)-high probability. Since a > ^/k/K, we obtain from (6.46), 

(^^.)-«+i-l^<C^«-^. 

The second term on the right side of (6.48), can be bounded by using (6.47). The first term on the right side 
of (6.48) is estimated by using the second inequality in (6.46). Thus, for any E satisfying (6.44), z G 2?l, and 
any A € X>Ao , we obtain that 



with (^, !^)-high probability. Thus 



Immiz) < Immfciz) + A(^) < , (6.49) 

Nrj 



with (^, i/)-high probability, for such E. By the spectral decomposition of H, we have 



lmm{z) = ^Yl 



N 

V 



N ^^if^a-Er + r,^' 
and we conclude that 

ImTO(z)>-^, (6.50) 

for some c > 0, if there is an eigenvalue in the interval [E — rj, E + rj]. Thus (6.49), implies, for any E, 
satisfying (6.44), that there is no eigenvalue in the interval [E — r],E + t]], with (^, i/)-high probability. 
To cover energies E > Eq, we use the following result: For a Wigner matrix W, satisfying 2.1, we have 

||W^||<2+^, (6.51) 

with (^, i/)-high probability. We refer, e.g., to Lemma 4.3. in [12]. Thus spectral perturbation theory implies 

\\H\\ < \\XV\\ + \\W\\ < 2 + + A, with {^,u)-high probability, covering the regime E > Eq. This concludes 
the proof. 

□ 
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6.4.2. Integrated density of states. In this subsection, we prove Theorem 2.18. Given the results on n(£^i,£'2) in 
Theorem 2.15 and the estimate on \\H\\ this is straightforward: 

Proof of Theorem 2.18. We assume that E is such that \E — Lx\ < \E — L2\. The other case is dealt with in the 
same way. Set 

with some ci large enough, such that UfdEi) = and n{Ei) = with i^)-high probability; see Lemma 6.7. 
Next, choose E > Ei, then from (2.23), we get, setting E2 = E and bounding £ < E — Ei + {ipn)^^N~^, 

\n{E) - nME)\ < (^^)^« (^1 + E - E, + M^^N-^^ . 

with i/)-high probability. Using our assumption on E and (6.52), we get 

i„(i^)-Mi^)i (1 + 1; + ^ + ^ 

with (^, :/)-high probability, for some C2 large enough. This estimate holds for any E and A. Uniformity is 
obtained with a lattice argument, we omit the details. □ 

6.4.3. Rigidity of eigenvalues. In this subsection, we prove Theorem 2.19. Recall the definition of the classical 
location ja of the a^^'-eigenvalue in (2.26). 

Lemma 6.8. There exists a constant C, such that, for all A e Paqj the following statements hold with (^, u)-high 

probability for some large enough c: 



i. //max{K^^,K^„} < {ipn)''^{-^ + j^), then 



|Ma-7a| < iVNf^ ( ^ ^ 



a. //max{K^„,K^^} > {ipn)'^{-^ + ^^573), then 



V]V SV3Ar2/3 7Vl/3a2/3 J ' 



where a :— minjo;, N — a}. 



Proof. We will focus on the eigenvalues /Ui, . . . ,/xjv/2- The other eigenvalues can be treated in a similar way. 
Define an event S as the intersection of the events on which the estimates 

\\H\\ < me^{\L,\,L2} + {mf°^ (-^ + ^) , (6.53) 

(see Lemma 6.7), and 

HE) - nf.{E)\ < i^.fo, (1 + ^ + 



N 

(see Theorem 2.18), hold, for any A < Aq and \E\ < Eq. We note that on S, we have ij.n/2 < K, for some K < L^. 
Let Ci > Cq. We use the dyadic decomposition 

21oe: AT 



{i,...,Ar/2}= y Uk, 



k=0 
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where 

Uo := {a < N/2 : + max{/7,„, 7„} < 2(^iv)^^«AiV-i/2} , 

Uk := {a < N/2 : 2''X{^Nf^^N-^'^ < + maxj/Xa, 7a} < 2*'+\^Nf-^\N-^/''} , (A; > 1) . 
By the definition of Uq and (6.53), we have 

on S, for a € Uq. 

For fc > 1, we find on S that 

^ = n/c(7a) = n(Ma) = n/c(/xa) + (¥'iv)^««0 + ]^ + ^) • (6.54) 
On S, and for a e Uk, we can bound the second term on the right side of the above equation as 



Af3/4 



where we used < \Li \ + /Xq,- Furthermore, we have on ^, for a G ?7fe, 

n/e(7a) + n/c(/^a) > c23'=/2 (^^)3C.€/2^3/2^-3/4 ^ 

where we used nfc{Li + a;) a;^/^, for < a; < 1^11+ ii'- Thus 

{<PNf°^0 (J^ + j^ + « n/e(7a) + n/c(/Xa) , 

which implies by (6.54) that 

UfM = n/c(7a) (l + O ((¥'Jv)-(^^-^°)«)) , 

on 5, for a € Uk- Using that n'f^{x) <^ {nfc{x)y/^ (|Li| +a;)^/^, for Li < a; < JC, we have +7a ^ +/ic«- 
Hence 

n'fcix) ~ n/c(7c<) , 

for any x between ja and ^a- Recalling that the density Pfc{x) is continuous, we conclude that, on S, for a G Uk, 

,\nfc{lJ'a) - n/c(7a)| 



IfJ-a -7a\ < C- 



"/c(7a) 



[a/Nyi^ \N 7V3/4 



where we used < + — 7o;| and /t-y^ ~ {a/N)'^/^. Next, since a = Nnfc{la) ^ -^(l-^il + 7a)^^^, we 
find for a e f/fe, {k > 1), 

3/2 



hence a ^ ^ 1/12 Using Young's inequality, we can absorb the last term on the right side of (6.55) into 
the left side and we obtain 

A2 a 

+ 



|Ma-7a| < (<PJV)^« ( ^1/3^^2/3 + 



a2/3iVi/3 ^ 



on S, for a G Uk, some C sufficiently large. The proof is completed by noticing that the event S has (^, z^)-high 
probability. □ 
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We conclude this section with the proof of Theorem 2.19. 

Proof of Theorem 2.19. We restrict the discussion to eigenvalues with a < N/2, the other eigenvalues are dealt 
with in the same way. From a/N = Ufd'^a) ~ (|Li| + 70)"^/^, we find that 

a< (¥'jv)^^(1 + A3/27V1/4)^ (6.56) 

if a is as in item i of Lemma 6.8. Combing the conclusions of items i and ii of Lemma 6.8 with (6.56) completes 
the proof of the theorem. □ 



A. Appendix: Free Convolution Measure and Stability Bounds 



A.l. Introduction. In this appendix, we discuss some properties of the free convolution measure, /U/c, defined 
through the functional equation 

mf.iz) = f ^^^'^'-^''^ , z = E + in^C+; (A.l) 
J-i Xv- z-mfc{z) 

such that Immfc{z) > 0, for 77 > 0; c.f., Equation (2.6). Here A > and we assume that fi is an absolutely 
continuous measure, with bounded and continuous density ^{v) such that supp/i = [—1,1]. For simplicity, we 
always assume that /x is centered, although this is not essential for our argument. 

To see that Equation (A.l) has a unique solution such that ImnifdE + irj) > 0, for ?7 > 0, one can choose 
r) > 2 first. Then it is straightforward to check that the right side of (A.l) is a contraction (in the sup- norm on 
the set of analytic function on the upper half plane with positive imaginary part). The fixed point equation (A.l) 
thus has a unique solution for 77 > 2. By analytic continuation, the solution extends to the whole upper half 
plane. We leave the details aside and refer, e.g., to [34]. 

A deep study of the equation (A.l), with slightly diff'erent conventions, can be found in [5]. One important 
result of [5] is the following: The measure /u/c is absolutely continuous with respect to Lebesgue measure, in 
particular, we have Trfifc{E) = lim^N^o Iin"^/c(^^ + i??)- For general probability measures (of bounded support), 
the support of iJ.fc may consist of several disjoint intervals, however, under our assumptions, the support of the 
fife is a single interval, say on [Li, L2], with Li < < L2; see Lemma A.l. We refer to [5] for a discussion of the 
general case. 

We are mainly interested in the behaviour of HfdE) and lmmfc{E+ir]), for S e M close to Li, L2 respectively. 
We distinguish the cases A < 1 and A > 1: 

For the former case, it was already pointed out in [5] (see also [31, 36]) that fife has a square root behaviour 
near L2, i.e., /U/c(L2 — k) ~ ^/k, k>Q, and similar for Li. 

For A > 1, we will restrict our attention to Jacobi measures, a special class of measures whose densities are of 
the form 

f,{v) = Z-\l + vni - vMv)x[-i,i]{v) , 

where a, b > —1, d € C^([— 1,1]) with d{v) > 0, v £ [—1,1] and the normalization constant Z is appropriately 
chosen so that becomes a probability density. Again, for simplicity, we will always assume that/U is centered. 
Note that we also admit exponents a, b smaller than zero, thus fi{v) — >■ 00 as w ^ ±1 is allowed. As it turns 
out, the square root behaviour at the endpoint of the support persists for A > 1, in case we have —1 < a, b < 1, 
respectively. However, if a, b > 1, there exists Aq > 1, such that for any A > Aq, we have /x/c(ii — k) ~ k^. For 
a precise statement; see Lemma A. 5. 

A. 2. Case A < 1. In this subsection, we choose A < 1. Adopting the proof of Proposition 2 in [36], we have the 
following result: 

Lemma A.l. Let fi be a centered probability measure supported on [—1,1]. Assume that /i has a continuous, 
strictly positive, bounded density ji{v) on (—1,1). Suppose that < A < 1. Then, there exits Li,L2 € M, with 
Li < < L2, such that the free convolution of fi with the semicircle law, fife, satisfies 

supp/U/c = [Li,L2] . 
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Moreover, denoting by ke the distance to the endpoints oj the support of Hfc, i.e., 

KE :=min{|i;-Li|,|E-L2|}, 



we have 



for some constant C > 1. 



C-^./^<^^fc{E)<C^, E&[LuL2], (A.2) 



We briefly outline how the proof in [36] can be adopted to our setting: We denote by mfc{z), z G C+, the 
Stieltjes transform of the free convolution measure /x/c- Define t := z + mfc{z) and consider instead of (A.l) the 
equation F{r) = z, where 

Fir):=r- [" M± , . e C+ . (A.3) 



I Xv — 

Note that limy\^olniF{x + iy) = —ir^{x) < 0, for x G (—A, A), since we have assumed that the density of fj. is 
bounded and continuous, and strictly positive in the interval (—1,1). Thus F extends to a function on M, which 
is continuous and bounded, except possibly at the point {±A}. As shown in [36], the end points, (Li), of the 
support of /x/c are characterized as the real valued solutions, r,, with \Ti\ > A, of the equation F'{t) = {Li are 
then obtained by solving Tj = Li + m/c(I/i)). Setting 

a point E is an endpoint of the support of /x/c, if H{r) = 1, \r\ > X, r = E + mfc{E) G R. Since A < 1 and /x is 
centered, we have from Jensen's inequality 

HiX) = [' > ^- = ^ > 1. (A.5) 

y_i(At;-A)2 A2 (jd/x(^^)(„_i))2 A2 

Here, the flrst inequality is strict since /x is absolutely continuous. Since H{t) is monotone decreasing (on R) as 

|t| — > oo, wc conclude that there are only two solutions ti, T2, such that the endpoints of the support of jifc, Li 
and L2 satisfy Li < —2 and L2 > 2. The square root behaviour of Hfc at Li, i.e., (A.2), follows as in [35]: It 
suflices to observe that F"{Ti) 7^ 0, thus by the inverse function theorem, we have, for 2; e C in a neighborhood 
of Li, F'^~'^\z) = Li + Ci^ z — Li{l + Ai{z — Li)) (such that Imi^~-'^(z) > 0, for z € C"*"), for real constants q 7^ 
and analytic functions Ai, with \Ai\ < 1 in a neighborhood of zero. This concludes our discussion on the proof 
of Lemma A.l. 

As an important corollary of the proof of Lemma A.l, we have the following stability bound already pointed 

out in [36]: 

Corollary A.2. Under the assumptions of Lemma A.l there exist constants C,c> 0, such that 

c < \Xv — z — mfc{z)\ < C , z = E + iri, 
for any Xv e (—A, A) and \E\ < Eq, < 77 < 3. 

Proof. For the upper bound, wc note that |m/c(-2^)| 1, as follows from considering the imaginary part of mfc 
in (A.l). For the lower bound, note that in a neighborhood of L,, Re(2;+m/c(2)) = 'ReT{z) = ReTi+0{\z—Li\^/^). 
Since \Ti\ > 1, \ Re(z + mfc{z))\ > 1, for \z — Li\ < e for a sufficiently small e > 0. For | Rez| > \Li\ + (e/2), 
the estimate is trivial. In the region not covered by the two preceding estimates, we must have Imr > c, thus 
Immfc + rj > c. The claim follows. □ 

Remark A.3. The idea in the proof of Lemma A.l can be applied to prove the analyticity of the density pfc 
of the free convolution measure /k/c. The analyticity of Pfc{z) is equivalent to the analyticity of mfc{z), or 
t{z) — z + mfc{z) for z G (Li,L2)- Since t{z) = F~^{z), it suffices to show that F'{t) ^ 0, for r e C+, where 
Im i^(r) = 0. 
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A. 3. Case A > 1. In this subsection, we choose for simplicity /z as a Jacobi measure, i.e., is described in terms 
of its density 

fiiv) = Z-\l + vni - v)''d{v)x[-i,i]{v) , (A.6) 

where a, b > —1, d e C^{[—a,b]) such that d{v) > 0, v £ [—a,b\ and Z is an appropriately chosen normalization 
constant such that /x is a probability measure. Below, we will assume, for simplicity of the arguments, that /x is 
centered, but this condition can easily be relaxed. 

Lemma A. 4. Let ^ be a centered Jacobi measure. Suppose that A > 1. //— 1 < a, b < 1, the results in Lemma 
A.l and Corollary A. 2 hold true. 

Proof. We can apply the same argument as in the proof of Lemma A.l. The only thing we need to prove is that 

for any sufficiently small e > 0, and a similar estimate; for H{—\ — e). 

From the assumptions, we find that there exist constants C, Co > such that ^,{v) > C(l — v)"" > Co(l — v) 
for any v G (0, 1). Let n := e^+^ '^o and choose e < 1/n. Then, we have 

rr,. -.^^ {l-v)dv Co r"^^t-{e/X), Co,, , 1, , 

Ji-(n-i)e/\ (Aw - A - e)2 A2 J^/^ t^ A^ n 

From the continuity of H, we can get the desired results. The same argument applies to H{—\ — e). □ 

For a, b > 1, we have the following result: 
Lemma A. 5. Let jj. be a centered Jacobi measure with a, b > 1. Define 

fi{v)dv \ ^ fj.{v)dv 



-1 ^-v 

Then, there exist Li <0 < L2 such that the support of Hfc is [Li,L2]. Moreover, 
i. if X < X2, then for < k < L2, 

C-^^< Hfc{L2- k) <C^, (A.7) 

for some C >1. 
ii. if X> X2, then L2 = X + (T2/A) and, for < k < L2, 

C-'^k}' < M/c(i2 -k)< Ck^ , (A.8) 

for some C > 1. 
Moreover, for < E < Eq, < r] < 2, z = E + ir], v e [-1, 1], 

\Xv — z — mfc{z)\ 

remains bounded from below in case ii uniformly in z and v, but in case i, it can be arbitrarily small as v ^ 1, 
E = L2, and 77 — >■ 0. 

Similar statements hold for the lower endpoint Li of the support of fife, with T2 and X2 replaced by 
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Proof. We first note that < A2,T2 < cxd for b > 1. Since fx{v) > 0, for v G (—1, 1), /i/c is supported on a single 
interval. Consider now 

When A < A2, we may follow the proof of Lemma A.l to prove the claims in i. 

We now choose A > A2. Then, -ff(A) < 1 and the curve on which ImF(r) = does not connect with the real 
axis on A, A]. Since this curve cannot end at some point where F is analytic, we can conclude that the curve 
approaches to the spectral edges, —A and A. When r = A, we have 



z = T- mfc{z) = X~ j 

To prove (A.8), let 



^ IJ.{v)dv _ . T2 

A^-^ + T 



T2 

T = X — Xk + iXy , z = X + — — k + ir] . 

A 



1^ uMdv y n(v)dv . . , „^ 

Ay-r7 = Imm/,(z) = Im / = ^ \ \ . (A.IO) 



Considering the imaginary part of m/c, we obtain 

^' - , 

Xv-T X J_i {v - 1 + k)"^ +y 
We claim that in the limit ?? \ 0, 

+ (A.ll) 

for K,y <^1. 

For the upper bound, we consider first the case y < k: Let e = min{l/2, (A^/A2) — 1}, then we have 

^J.^{v-l + kr+y^-y[J_, + A_s,-.fc + + J {l-v-k^+y^- ^ ■ ■> 

The first term in (A. 12) can be estimated as 

<(i + ^)a^.. 

Here, we used that v <1 — 8e~^k implies that 1 — v — k > (1 — e/8)(l — v), hence 



(l-w-fc)2-(^ s) (1-v)^ - 2) (1-v) 



(l_-y_fc)2-v 87 (l-t;)2-V 2J {1-v)'^ 
The second term in (A. 12) can be estimated as 

'■^-''-y a(v)dv ^ r^-''-y k^dv 



A-8e-. {1-v-kr+y^-^y L,,-., ii-v-ky-^" ■ 



The third term in (A. 12) can be estimated as 

, , f^^', ^ O, n- &±4± ^ Ciy . t)' 

Ji-k-y {v-l + ky + y^ Ji-k-y r 

The last term in (A. 12) can be estimated as 



y 



[' n{v)dv {k-y)^dv ^ f'^ {k-yfdw ^ 
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Thus, as \ 0, we have that 
Since 

^~XV^2) ^2-TVA|"^"2j - 2A ' 

we obtain 

y<C{k + y)'^, 

provided y < k. 

When y > k, we decompose the integral in (A. 12) as 



(1 - ?; - + 2/2 • 

The first term is again estimated as in (A. 13). The second term can be estimated as 

j\-s,€-^k K'- ~ ^ ~ +y Ji-se-ifc y 

Following the argument we used for the case y < k, we find the relation y < Cy° in this case. For sufficiently 
small y, this is impossible, so this case does not happen. 

To complete the proof of (A. 11), we need a lower bound: Observe that 

yj.Av-i + ky + y'-''l.-k-y{v-i + kY + y--^y y- ^""-'^^ 

and (A. 11) follows from (A. 10) and that k > y. When y,k <^ 1, (A. 11) implies that k ^ y and since y Cjifc 
as 7/ 0, we have (/x/c) ^ y ^ k^. 

To compare k and k, we consider the real part of m/c and get 

^ li(v)dv 1 (v — 1 + k)ii(v)dv 



K-Xk-^=RemfJz)=Re , ' ^ - - 

A ' J-iXv-T A y_i {v - 1 



+ fc)2 + 2/2 • 



From the definition of T2, we find that 

_ 1 nv-l + k)iJ,{v)dv iJ,{v)dv\ _ 1 /■^ iJ,{v)dv - 1) + P + j/2 

~xj-l[{v-l + k)^+y^ ^T^j ~ Ay_iT^' (^;-l + fc)2+t/2 

We now separate the integral and estimate each term as in (A. 14) and (A. 12). We then get 



kf^{v)dv ..k^k + yt 



2 



_1 {V -1 + fc)2 + y2 y 

and 

1 li{v)dv e + y^ e + y 

Xj_il-v {v-l + ky + y^^ y ^ ■ 
Recalling that y ~ k^, when y, <C 1, we find that k — Xk — 0{k). Therefore we get 

fifc{L2 - k) y k"" K° , 

as /v \ 0. Finally, it is easy to see that \Xv — z — m/c(2;)| is not bounded from below: Choosing z = L2, we have 
lm(rn/c(L2)) = 0, but Re(Aii — L2 — mfc{L2)) = Xv — X. This proves the claims in ii. 

□ 
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A. 4. Square root behaviour of m/c and further stability bounds. In this subsection, we prove that the 
Stieltjes transform m/c inherits the square root behavior from fife- 

Lemma A. 6. Assume that /x/c has support [Li,L2] and satisfies 

C-^^/k< lifc{L2- k) <C^, (A.15) 

0< K< L2, C >1. Then, 

i. for z = L2 — K + iT], with < k < L2 and < ry < 2, we have, C > 1, 

C~^^K + r\ < lmmfc{z) < C\/k + t] . 

ii. for z = L2 + K + ir], with < k <1 and < r] <2, we have, C > 1, 

The analogous statements hold for z = Li± K + 'vq. 
Proof. We start with the claim i: Notice that 

>l tl/'/c(-i-) 



{x - L2 + kY + rf 



To prove the lower bound, consider the following cases: 

Case 1. When K,ri < 1/2, computing the integral from x = L2 — k — 2r] to x = L2 — k — ri, we find from (A.15) 
that 

Case 2. When n > 1/2, 77 < 1/2, wc obtain (A.15) that 



Case 3. When 77 > 1/2, we have a bound 



' (X — L2 — k) + J rj 



This proves the lower bound. To prove the upper bound, we consider the following cases: 
Case 1. When r] < k < 1/2, from (A.15) we have 



Imm/c(-2) = j dnfc{x) 



r] 



{X-L2 + k)2 + 7?2 



<Cr;/ ^^^^dy + C^^+lj + C?^ \dy<CV^^. 

J Tj y J rj y 

Case 2. When k < rj < 1/2, calculation similar to Case 1. proves the same bound. 
Case 3. When k> 1/2, we have 

Imm/e(^) <C [ , <C< G^f^. 

J [x- L2 + kY + r)^ 

Case 4- When r] > 1/2, we have a trivial bound 



lmmfc{z) < \mfc{z)\ < - < C^k + rj . 
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This completes the proof of statement i To prove ii, we proceed similarly: 

Case 1. When k> r], computing the integral from x = L2 — k to x = L2 — 2k, we get 



For the upper bound, we find that 

rjdfifcix) 



Immfc{z) = 



(x — L2 — k)^ + rj"^ 



-dx + Ct] 



-L, {X-L2f 



dx 



J K 



Case 2. When k< r], computing the integral from x = L2 — (??/2) to x = L2 — r], we obtain 

L2-ri 



J {x-L2-KY + rf Jl 

For the upper bound. 



L2-(r,/2) n 



^-4dx>cvv>^ 



Imnif, 



,(.) = / 



{x - L2- kY + rf 



L2-V 



dx + Cr] 



dx 



<C^ + Ct] f ^<cVv 



<cvv<^^. 

\/r] + K 



□ 



Finally, we show that |1 — ii2|/|-R3| ~ \/he + t]\ see (3.31) for the definition. For simplicity, we assume that 
]ifc has a square root behavior at both endpoints of its support. 

Lemma A. 7. Assume that Hfc has support [Li,L2] and satisfies 

C-^y/^<]ifciE) < C^, 

with C > 1, where ke '■= min{|£^ — Li|, \E — -L2I}, denotes the distance to the endpoints of the support of ]j,fc- 
Moreover, assume the stability bound 

c < \X - z - mfc{z)\ < Co , 

with Cf),c > 0, for \E\ < Eq and < ry < 3. Then, there exists a constant C > 1 such that for any \E\ < Eq, 
<r]<3, 



C-'^s/KE + r]< 



d]j,{v) 



{\v- z- mfc{z)Y 



d]2{v) 



{\v- z- mfc{z)y 



< C^/ke + r] . 



Proof. Since c < \\v — z — to/c(z)| < C, it suffices to show that 



1 



d]i{v) 



{\v- z- mjc{z)Y 



< C^/ke + r]- 



We first consider the following decomposition: 

d]j,{v) 



1 



< 



1 



(Aw — z — mfc{z)Y 
d]j,{v) 



/d]i{v) f dn{v) f dfxi 

\Xv - z - mfc{z)\'^ J \Xv - z - mfc{z)\'^ J {Xv - z - 



dfi{v) 



mfc{z)Y 



53 



For the first term in tlie riglit side of the decomposition (A. 4), we have 



Imnifciz) 



< 



Im.{z + mfc{z)) lia{z + mfc{z)) C^/ke + rj 



a E € [Li,L2] and 



iiE€ [Li,L2Y. Since 



Im mfc{z) 



< 



lm.{z + mfc{z)) Im(z + m/c(z)) Crj/ ^ke + ?7 



< Cs/ke + ti, 



I Re(At; — z — m/c(2:))|,Im(2; + mfc{z)) <\z + m}c{z)\ + A < C , 



we also find that 



/ 



dfi{v) 



= 2 



\Xv ~ z — mfc{z)\'^ J {Xv — z — mfc{z))^ 

(Im(z + mfc{z)))^ + iRe(Av ~ z — nifdz)) ■ lm{z + mfc{z)) 



I 



\Xv — z — mfc{z)\^ 
< C lm{z + mfc{z)) < C\/ke + V ■ 



dfi{v) 



Thus, 



I {Xv 



diJ,{v) 



z - mfc{z)f 



< Cs/ke + V , 



which proves the upper bound. 

For the lower bound, we first consider the case E e [Li, L2]: Observe that 



Im 



Hence, 



{Xv - z- mfc{z)y 



djiiv) 



I Rc(Ar - : - ■ Iiu(-- + nij^jz^) 

\Xv- z- mfc{z)\'^ 



dii{v) 



> C lm{z + mfc{z)) > Cy/KE + f] . 



j {Xv 



z-mfc{z)y 



> 









diJ,{v) 



{Xv - z - mfc{z)Y 



> C^/ke + ti, 



for E & [ii, L2\. In case E e [Li, ^2]'^, we obtain a lower bound from 



J {Xv-z- 



dii{v) 



and 



1 



-/ 



dii{v) 



= 1 - 



mfc{z)Y 
lm.mfc{z) 



> 



I 



dn{v) 



\Xv - z- mfc{z)\'^ 



(A.16) 



> 



\Xv - z - mfc{z)\'^ lm{z + mfc{z)) lm{z + mfc{z)) r]/\/>^E + V 

This completes the proof. 



= Cy/KE +??. 



□ 
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